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[1] We have constructed new apparent polar wander paths (APWPs) for major plates over 
the last 200 Myr. Updated kinematic models and selected paleomagnetic data allowed 
us to construct a master APWP. A persistent quadrupole moment on the order of 396 of 
the dipole over the last 200 Myr is suggested. Paleomagnetic and hot spot APW are 
compared, and a new determination of “true polar wander” (TPW) is derived. Under the 
hypothesis of fixed Atlantic and Indian hot spots, we confirm that TPW is episodic, 
with periods of (quasi) standstill alternating with periods of faster TPW (in the Cretaceous). 
The typical duration of these periods is on the order of a few tens of millions of years with 
wander rates during fast tracks on the order of 30 to 50 km/Myr. A total TPW of some 30° 
is suggested for the last 200 Myr. We find no convincing evidence for episodes of 
superfast TPW such as proposed recently by a number of authors. Comparison over the 
last 130 Myr of TPW deduced from hot spot tracks and paleomagnetic data in the 
Indo-Atlantic hemisphere with an independent determination for the Pacific plate supports 
the idea that, to first order, TPW is a truly global feature of Earth dynamics. Comparison 
with numerical modeling estimates of TPW shows that all current models still fail to some 
extent to account for the observed values of TPW velocity and for the succession of 
standstills and tracks which is observed. ^ INDEX TERMS: 1527 Geomagnetism and 
Paleomagnetism: Paleomagnetism applied to geologic processes; 3040 Marine Geology and Geophysics: Plate 


tectonics (8150, 8155, 8157, 8158); 8120 Tectonophysics: Dynamics of lithosphere and mantle—general; 
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1. Introduction 


[2] Analysis of the fossil magnetization preserved in 
rocks is the basis for constraining such diverse geophysical 
problems as dynamo generation in the Earth's core, plate 
kinematics and paleogeographic reconstructions, and mantle 
dynamics leading to true polar wander (TPW). The second 
half of the twentieth century saw the advent and consol- 
idation of plate tectonics: paleomagnetic measurements on 
lava and sediments coming mostly from continental areas or 
from hot spot volcanics demonstrated continental drift and 
could be blended into apparent polar wander paths (APWP); 
oceanographic exploration led to the discovery of seafloor 
spreading related magnetic anomalies and transform faults, 
allowing the construction of kinematic models for each 
individual ocean basin. APWPs and kinematic models were 
developed from the '60s to the '80s under the key assump- 
tion that, when averaged over a sufficient amount of time, in 
excess of a few thousand years, the Earth's magnetic field 
could be described accurately by an axial centered dipole. 
Of course, both types of largely ocean- and continent-based 
data could not be independent. Over a decade ago, we 
proposed [Besse, 1986; Besse and Courtillot, 1988, 1991, 
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hereinafter referred to as BC91] to blend the two approaches 
and to use then available paleomagnetic data from North 
America (NAM), Africa (AFR), Eurasia (EUR) and India 
(IND) and kinematic models from the Indian, central 
Atlantic and North Atlantic oceans into a single “synthetic” 
APWP that could next be transferred to any desired plate. 
On the basis of a selection of 111 poles, the paths were 
defined in 20 Myr windows extending back to 200 Ma. This 
first-order analysis provided a reasonable check of the 
consistency of individual plate paleomagnetic data and 
kinematic models with the rigid plates and geocentric dipole 
hypotheses. It also filled gaps in the available APWPs from 
individual plates. However, the study suffered from a num- 
ber of limitations in space and time resolution and was not 
truly global since a number of major plates (South America- 
SAM, Australia-AUS, Antarctica-ANT and the Pacific-PAC) 
were not included. 

[3] Slight departures from a purely centered dipole field 
had been noted as early as 1970, when Wilson [1970] argued 
for a far-sided and right-handed distribution of virtual geo- 
magnetic poles during the Cenozoic. A number of analyses 
confirmed that when averaged over the last few million years, 
an axial quadrupole component is detectable, with an ampli- 
tude on the order of 3 to 6% of the axial dipole [e.g., Coupland 
and Van der Voo, 1980; Merrill and McElhinny, 1983; 
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Constable and Parker, 1988; Schneider and Kent, 1990b; 
Quidelleur et al., 1994; Carlut and Courtillot, 1998; Johnson 
and Constable, 1997]. Livermore et al. [1983, 1984] tried to 
extract such a quadrupolar term in a worldwide paleomag- 
netic database going back 200 Ma. Besse and Courtillot 
[1991] argued that, for the period prior to 5 Ma going to 
200 Ma, a significant quadrupolar term could not be extracted 
unequivocally from the data available at that time. 

[4] On the other hand, comparison of paleomagnetic and 
hot spot APW led Besse and Courtillot [1991] to infer that 
significant true polar wander, amounting to more than 20°, 
had occurred in an episodic, irregular way in the last 200 
Myr. Recent advances in seismic tomography and dynamic 
modeling have raised fascinating new questions on the 
existence and role of slab avalanches, plume events and 
the viscosity structure of the deep mantle [e.g., Machetel 
and Weber, 1991; Ricard et al., 1993; Tackley et al., 1994; 
Steinberger and O’Connell, 1997; Bunge et al., 1998; 
Evans, 1998; Richards et al., 1999; Greff-Leffiz, 2001; M. 
Greff-Lefftz and P. Bunge, personal communication, 2000] 
(see, e.g., review by Tackley [2000]). Reliable measures of 
TPW are required to constrain these models. 

[5] Paper BC91 was based on a database compiled and 
updated in our institute [Besse, 1986], prior to availability 
of the global paleomagnetic computerized database of 
McElhinny and Lock [1995]. In its June 1999 version, 
the GPMD includes all paleomagnetic data published up to 
the end of 1997. As will be seen below, and still retaining 
the selection criteria used in BC91, the total amount of 
reliable data has been increased by a factor of more than 2. 
Also, significant methodological progress, most notably 
regarding satellite altimetry data analysis, has led to much 
improved description of fracture zones and continent/ocean 
boundaries and hence better reconstructions of past plate 
motions. Improvement in normalization and correlation of 
timescales has resulted in data sets such as the Global 
Seafloor Anomaly Chart [Royer et al., 1992] in which the 
marine geophysics community has integrated most avail- 
able data into a single coherent synthesis. 

[6] As a consequence of publication of these large sets of 
new data, we have felt it useful to update our BC91 paper a 
decade after its publication to try and evaluate whether 
improved plate reconstructions, APWPs, constraints on 
field geometry and mantle dynamics could be obtained. In 
order to decrease the risk of generating artifacts due to 
limited site distribution [see, e.g., Quidelleur et al., 1994], 
we have updated the NAM, EUR, AFR and IND data sets 
and have added data from AUS, ANT, SAM and GRE 
(Greenland). We have also included some data from ODP/ 
DSDP drilling sites and two skewness poles coming from 
analyses of marine magnetic anomaly profiles from the 
Indian Ocean (potential problems when using this type of 
data are discussed below). This results in much improved 
geographical coverage both in latitude and longitude (Fig- 
ure 1). The case of the Pacific plate and data is left to the 
discussion section. 


2. Extended Database 
2.1. Paleomagnetic Data From Continents 


[7] Faster paleomagnetic data search and treatment is now 
made easy due to the availability of the computerized Global 
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Paleomagnetic Database (GPMD) of McElhinny and Lock 
[1995]. We have used the June 1999 version (GPMD V3.3). 
In the process of using and searching the database, we have 
uncovered a few errors and inconsistencies. The more 
frequent occurrences were older data later re-used in a larger 
study or an update without any notice of this redundancy in 
the database. In order to avoid confusion in the database, and 
to allow others to easily check our computations, we have 
preferred not to alter the database version we started with, 
but rather send a list of the problems to the managers of the 
database for future revisions, equally accessible to all. We 
have, however, checked that these errors and inconsistencies 
were sufficiently few and small that their integration in the 
master APWPs resulted in changes of mean pole positions 
always less than 1°, or in mean age less than 1 Myr. As a 
result, these changes would hardly alter the Figures. 

[8] The key question is which selection criteria are 
applied to sift the base and extract a robust high quality 
subset that allows one to address a given problem. As was 
already the case in BC91, we have decided to retain data 
from both magmatic and sedimentary rocks. In a similar 
recent study of TPW, Prévot et al. [2000] chose to retain 
only magmatic rocks, which they believe to carry a more 
reliable magnetization, namely a thermoremanent magnet- 
ization (TRM). The question of the relative reliabilities of 
sedimentary and volcanic magnetizations 1s an old one in 
paleomagnetism. For instance, deepwater sediments can 
have significant compaction errors [e.g., Celaya and Clem- 
ent, 1988; Tarduno, 1990] but there are none in our data set. 
We believe that there can be as many problems with 
volcanic rocks as with sediments (e.g., remagnetizations 
versus inclination shallowing, too short versus too long 
integration time for recording of magnetization) and in any 
case it is interesting to compare results obtained in either 
case, as will be done below. 

[9] For the sake of comparison, and because we believe 
the end result to be sufficiently reliable, we have retained a 
selection procedure similar to BC91. Applying selection 
criteria automatically, 1.е., writing selection filters, is quite 
straightforward with a properly computerized database. The 
BC91 criteria were the following: (1) at least 6 sites and 36 
samples per study; (2) a 95% confidence interval less than 
10? in the Cenozoic and 15? in the Mesozoic; (3) evidence 
for successful alternating field and/or thermal demagnet- 
ization (i.e., Demagnetization Code equal to or larger than 2 
in the McElhinny and Lock [1995] terminology); (4) dating 
uncertainties less than 15 Myr; and (5) absence of remagne- 
tization ensured by the fact that differences between mag- 
netization and stratigraphic ages should be less than 5 Myr 
and by rejecting poles with a negative fold test or negative 
reversal test. 

[10] These simple criteria are of course to a large degree 
arbitrary and are based on our experience in trying to 
eliminate most problematic studies. This is discussed in 
BC91 and compared to different selection criteria used by 
other authors. Part of the justification comes a posteriori 
from a comparative examination of the results. 

[11] One of the key aspects of our selection process has 
been an attempt to identify mobile zones from which data 
should be excluded from what should in principle be rigid 
plates. This is important not only for plate reconstructions, 
but also for field geometry analysis: Carlut and Courtillot 
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Figure 1. (a) Master apparent polar wander path for Africa from the Present back to 200 Ma, with 
associated ellipses of 95% confidence in shaded gray (averages every 10 Myr, with a 20 Myr sliding 
window); actual mean ages for each time window are indicated. (b) Distribution of individual 
paleomagnetic poles and sampling sites; dots: all poles used in this study to construct the master APWP, 
transferred to South African coordinates according to their age; stars: site locations reconstructed to 
appropriate time in South African frame (solid star, Northern Hemisphere, open star, Southern 
Hemisphere). Histograms from the Present back to 200 Ma for (c) number of data, (d) confidence 
intervals at the 95% level, and (e) Fisher's concentration parameter K. 


asuoor] suoururo;) oAnvar;) o[qeorjdde әці Aq рәшәло8 әле səpe WO ÉN Jo soni 10] A1e1qr IUHUO Кә uo (suonipuoo-pue-suroj/uio» Kop Kre1qr]ouruo//:sdiu) suonipuo;) pue SUHOL eq ess '[TZOZ/T L/Fc] uo AreaqrT ouuo Кә 000008 f000£/6c0T 0 L/I0p/uto2"o[TA" AxeAqrpoupuosqndn8//:sdnq woy popeopu«oq * TT ‘TOOT 02951 


EPM 6-4 


[1998] noted, for instance, in a study of the mean geo- 
magnetic field over the last 5 Myr that inclusion of tectoni- 
cally rotated Ethiopian data produced a very erroneous 
mean field model. Some tectonic areas, such as mountain 
ranges (Himalayas, Alpine ranges) or rift zones (Afar in 
Africa), are of course rather obvious and sites from these 
regions were eliminated. Others zones, such as the Colorado 
plateau in North America or the African rift south of the 
Afar depression are less straightforward. And actually, rocks 
often outcrop simply because of some amount of tectonic 
deformation. Our task is then to test whether the extent of 
deformation requires that the data be rejected or not. A 
major source of error is linked to local tectonic rotations 
about vertical axes: crustal deformation at a large scale or at 
depth may lead to small block rotations and many paleo- 
magnetic studies fail to take into account the plunge of fold 
axes (the resulting error may exceed 15°). Detailed analysis 
of individual poles and comparison with synthesized 
APWPs may allow one to spot remaining outliers with a 
previously undetected tectonic origin. In some cases, we 
have used another method: some poles that where suspected 
to have undergone rotations without significant poleward 
motion have been integrated using a method derived from 
McFadden and McElhinny [1988]. We treated these poles 
using only inclination data, 1.e., as small circle constraints, 
as we also did for data coming from DSDP Legs (see 
below). These poles (11 studies in western North America, 
one in South America and one in Europe) are marked in 
Table 1. 

[12] Another major problem is to determine to what 
extent large plates, such as Eurasia or Africa, can be 
considered as rigid. The ocean-based kinematic syntheses 
of Müller et al. [1993] imply that Africa and South America 
need to be subdivided into a small number of rigid subplates 
that have undergone some (slight) amount of differential 
rotation since the time of breakup. Recent surveys of South 
Asian data [Cogné et al., 1999; Yang and Besse, 2001] have 
emphasized differences between the Cretaceous segment of 
the Chinese APWP and the synthesized APWP of BC91 for 
Eurasia at the same period. Cogné et al. [1999] concluded 
that "rigid" Eurasia may have actually undergone some 
amount of internal deformation, possibly due to the relative 
motions of three rigid subplates. Therefore, and contrary to 
our analysis in BC91 where China was assumed to have 
been rigidly attached to Eurasia (hence Europe) since the 
Cretaceous, we have excluded Chinese data from the 
present compilation. This final data set comprises 221 poles 
(Tables la—1g). 


2.2. Oceanic Data 
2.2.1. Data From Skewness of Marine Magnetic 
Anomalies 

[13] The skewness of a (usually marine) magnetic anom- 
aly is the phase shift required to restore the shape of the 
measured anomaly to the one which would have been 
observed, had the anomaly been created at the geomagnetic 
pole [Schouten and McCamy, 1972]. The anomaly is then 
“undistorted”, i.e., has the same symmetry elements as the 
magnetic source body. This is particularly useful for accu- 
rate determination of the location of the center or edges of 
individual crustal blocks corresponding to a given chron. 
The value of skewness depends on the inclination. and 
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declination of the original remanence of the block, on the 
present field direction (inclination and declination) and on 
the azimuth of the magnetized body (1.е., generally of the 
ridge where the crustal block was generated). A great circle 
of possible pole positions is deduced from the skewness 
value in the hypothesis of a dipolar centered geomagnetic 
field. The intersections of a set of such semicircles derived 
from the same anomaly (chron) at distant locations on the 
same plate in principle allow the determination of a virtual 
geomagnetic pole [e.g., Schouten and Cande, 1976; Gordon 
and Cox, 1980]. 

[14] The analysis of the skewness of marine magnetic 
anomalies generated by seafloor spreading is therefore a 
way of obtaining pole positions for oceanic plates, for 
which other kinds of data may be lacking. The method 
has been extensively applied by R. Gordon and colleagues 
[e.g., Petronotis and Gordon, 1999] to the case of the 
Pacific plate. Because of problems in connecting that plate 
to the others through a plate circuit, we have not included it 
in the set used for construction of synthesized APWPs; we 
return to this question in the final section of this paper. We 
note, however, that the theoretical potential of the skewness 
method for the determination of VGPs from purely oceanic 
plates is unfortunately hampered by strong limitations due 
to the parasitic effect of neighboring sources and/or to the 
tectonic tilt of blocks generated close to the spreading 
center, leading to strong artifacts: Cande [1976] discusses 
the problems of “anomalous skewness”, which are further 
elaborated on by Petronotis et al. [1992] and Dyment et al. 
[1994]. 

[15] In the case of the Indian Ocean, careful analysis 
[Dyment et al., 1994] of conjugate marine magnetic anoma- 
lies from the Carlsberg ridge and the Wharton Basin, and 
the Central and South-East Indian ridges, leads to the 
conclusion that anomalous skewness becomes negligible 
when the spreading rate is faster than 50 km/Myr. The 
determination of skewness from anomaly 34 (84 Ma), 1.е., 
the time when India began its northward motion, up to 
anomaly 21 (47 Ma), when a sudden velocity drop resulted 
from the India-Asia collision [Patriat and Achache, 1984] is 
possible because the spreading rate is far above this critical 
value. Dyment et al. [1994] found that anomalies 21, 24 and 
29 were suitable for the determination of reasonably accu- 
rate paleolatitudes; they found good agreement between the 
great circle of possible pole positions for these anomalies 
and the BC91 APWP. 

[16] However, Dyment et al. [1994] did not publish such 
great circles for anomalies 25 and 26 which are prominent in 
the Indian Ocean. 7orcq [1997] and Е. Torcq and J. Besse 
(manuscript in preparation, 2002) filled this gap, by analyzing 
profiles from the Somali, Madagascar, Arabia, Crozet, and 
Central Indian basins. Anomalies 25 and 26 (61.5—57.2 Ma 
following the Harland et al. [1989] timescale) are easy to 
recognize and their shapes are well determined. Anticipat- 
ing results given later in this paper, the intersections of 
skewness derived circles for the Indian and African plates 
yield a pole in excellent agreement with data from other 
continents. There are unfortunately no other sets of usable 
skewness data from the plates involved in our synthesis. 
2.2.2. Oceanic Data From Boreholes 

[17] Paleolatitudes can also be inferred from inclination 
measurements performed on unoriented cores [e.g., Cox and 
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Year 


RNO 


Rock Type 


Tests 


Plong Dp 
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Age 


Rock Name 


Table 1g. (continued) 
North Carolina, northwest dikes 
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с D MO — 00 C1 CI CY оо — е ч © с 
RRRRRRRRRRE RF 98 RF Gordon, 1984]. Despite the huge number of cores which 
have been drilled in recent decades in the various steps of 
eae "P international drilling projects, only a few reliable data can 
neg$4ü9$59 K а За be used for the determination of paleolatitudes. Peirce 
[1976, 1978] provided early compilations of these data, 
А and selected the most reliable studies. One of the main 
$ $ problems encountered in oceanic cores is the absence of 
- К $ E E с E 2244 orientation: declination is not available, unless one uses the 
i 4 ЕЕЕ Е PES 8 8 present Earth's field to reorient the core. Inclination alone 
$3 S82 ge 8s BB $ allows only the determination of small circles of equal 
goes g "EP g ££ А А PE paleocolatitude on which the pole must lie. Another problem 
$ C a E E 5 E 2 B СЕЕ B is the underestimation of inclination due to arithmetic 
E E E В ZE В E БЕ Р EC B 33 D averaging (but see Enkin and Watson [1996]). In order to 
patio B overcome this problem, Peirce [1976] used an unpublished 
5 method of Cox, which also gives an estimate of the standard 
2253 wamme > б xw error (Sm). We have selected those data corresponding to 


о values of Sy, less than 5°. A minimum amount of 36 samples 
in a study was used, as was the case for terrestrial paleo- 
magnetic poles. 

[18] In addition, as regards data postdating the Peirce 
[1976, 1978] compilations, we have used the paleomagnetic 


benSeeer & F E A3 studies of oceanic cores drilled during campaigns of the 
Ocean Drilling (ODP) and Deep Sea Drilling (DSDP) 
L006508g8-9g - 9 шо Projects in the Indian (Wharton Basin and Ninety-East 


Ridge) and South Atlantic (South American plate) Oceans. 


Be ao MA E Eu x. A Aa We used the Leg 73 paleolatitudes published by Tauxe et al. 
Oo uw) — «= M0 pe . 
DE "DT ward E S en pcm n [1983], and the Leg 115 data of Schneider and Kent [1990a] 
and Vandamme and Courtillot [1990]. The latter studies give 
шуша a а 22 paleolatitudes for both the Indian and African. plates. All 
EA SS тк these results are listed in Table 2b. The distribution of these 
EE dv d © data, altogether representing 19 determinations from 0 to 
AACTCGAAANT m m ws 120 Ma, is also plotted in Figure 1а. 
Soe me Ea dv Ss © 2.3. Combining Dati Types o 
[19] The space and time distributions of selected paleo- 
aM aan = = ; | i 
S E 8 ЖЕШ E z z 3 а magnetic data are shown in map (Figures la and 1b) and 
= = = histogram (Figures 1c—1e) form. A total number of 242 
ооч атама w ^ ^9 independent data have finally been retained. This final data 
шс“? е з aS set comprises 221 "classical" poles, two poles determined 
using marine magnetic anomaly skewness data and 19 
ARAMSZwHSQRg Xi X ме “inclination only” data from oceanic cores. The total 
number of data is thus increased by more than 2 with 
58855555538 Б a 55 respect to the BC91 database, which comprised only 111 


poles. The map (Figure 1b) shows all poles (as dots) and the 
corresponding sites (as stars) in (South) African coordi- 
nates, with sites restored to their relative positions at the 
appropriate time. We see that geographical coverage is quite 
reasonable; it is much improved over BC91. The histogram 
(Figure 1с) shows that a large number of data (more than 
30) is available for the 0—30 and 50—80 Ma time windows. 
By contrast, the 140—150 Ma time window contains only 
eight studies. The mean number of data in each 10 Myr time 
window is on the order of 15. As a matter of comparison, 
the corresponding mean for the recent Prévot et al. [2000] 
study is on the order of 4, and there are less than 5 data in 
twelve 10 Myr time windows. 


3. Kinematic Models 


[20] Despite the significant improvement in the paleomag- 
netic database, the number of data is still rather low if one 
wants to construct the APWP of a single plate based on its 


“See Table 1a footnote for explanation. 


>See Table 1f. 


Vermont, New Hampshire, White Mountain Volcanics 

New Brunswick, Canada, Caraquet dike 

Nova Scotia, Canada, North Mountain Basalt 

Pennsylvania, Culpeper Basin intrusives and baked sediments 
New Jersey, Newark Basin Hettangian red beds 

New Jersey, Passaic Formation, Preakness Basalt 


South Carolina Piedmont, Diabase dikes 
South Carolina Piedmont, Diabase dikes 
Combined Pennsylvania diabase 
Connecticut Valley Volcanics 


North Carolina, north-south dikes 
Newark Group 


Newark Supergroup Volcanics 
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Table 2a. Pole Positions Deduced From Skewness of Marine 
Anomalies 24 and 25 in the Indian Ocean? 


Mean Age Anomaly CN) (CE) ND 95%СІ K 
57.2 25 76.8 224.7 17 6.9 30.7 
61.5 26 76.4 200.6 14 6.2 48.2 


“Definitions: \(°N), pole latitude; ф(°Е) pole longitude; ND, number of 
determinations; 95%CI, confidence interval; K, precision parameter. 


own paleomagnetic poles only. As in BC91, we therefore take 
advantage of the fact that the APWPs from individual plates 
cannot be independent one from the other and should be 
related through plate kinematic models. And indeed, accurate 
up-to-date models are available for the Indian, North Atlan- 
tic, central Atlantic, and South Atlantic Oceans, allowing one 
to relate the seven selected plates over the entire period of 
interest here with reasonable confidence. A possible motion 
between Eurasia and North America [Van der Voo, 1993] may 
constitute a source of error during Jurassic times. These 
models can be used to transfer all paleomagnetic data to a 
common, single reference frame, arbitrarily taken here to be 
the southern part of the African plate. 

[21] Müller et al. [1993] have published a new set of data, 
in which plate motions with respect to hot spots are 
computed based on combined Atlantic and Indian Ocean 
hot spot tracks between the Present and 130 Ma. The initial 
relative plate motion parameters are those quoted in the 
database of Royer et al. [1992], and are computed at any 
given age for the same anomaly. We have taken advantage 
of these compilations, both because of their unifying fea- 
tures, and because one of the aims of this paper is indeed to 
compute relative motions between the paleomagnetic and 
hot spot reference frames. Müller et al. [1993] extensively 
discuss the implications of their plate motion models. Two 
important features of their reconstructions are the internal 
deformation of the African and South American plates 
during breakup phases in the Cretaceous, and a new model 
of the final breakup of Gondwana, in which India rifted 
away from Antarctica at the same time it rifted away from 
Australia (1.е., between chrons M10 and M11). 

[22] On the basis of the Royer et al. [1992] and Müller et 
al. [1993] kinematics, we have recomputed the relative 


Table 2b. Selected Data From ODP and DSDP Legs? 
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motions at the time of each anomaly between 0 and 
130 Ma, by simply subtracting the absolute motion of 
(South) Africa with respect to hot spots from the motion 
of other plates (Table 3). The finite poles of rotation which 
are not given by Müller et al. [1993] are those of Müller and 
Roest [1992] for the motion between Europe and North 
America, and those of Srivastava and Tapscott [1986] for 
Greenland and North America. They were combined with 
the North America versus South Africa motion to obtain 
total motion with respect to South Africa. The Arabia to 
Africa rotation poles are those of Le Pichon and Gaulier 
[1988]. For periods earlier than 130 Ma, the South America 
to South Africa kinematic parameters are those of Nürnberg 
and Müller [1991], the Australia to East-Antarctica to South 
Africa rotation poles are those of Royer and Sandwell 
[1989], using respectively a recomputed fit and a fit after 
Lawver and Scotese [1987]. Prior to 130 Ma, the Mada- 
gascar to South Africa motion uses the parameters from the 
Global Isochron Chart of Royer et al. [1992], with a final fit 
from Lawver and Scotese [1987]. 

[23] For each plate, paleomagnetic or skewness poles and 
small circles based on DSDP/ODP inclination data were 
transferred onto (South) Africa. For this, we interpolated a 
Eulerian pole of rotation between the two published finite 
rotation poles with ages bracketing the estimated age 
respectively of the pole and drill hole or crustal block 
corresponding to the individual datum. Errors incurred in 
such reconstructions are unlikely to exceed 2? (see Molnar 
and Stock [1987], discussion by Besse and Courtillot [1988, 
1991], and Acton and Gordon [1994]). We must be cautious 
in case the timescales used by different authors are not 
identical and therefore require some discussion. McElhinny 
and Lock [1995] used the Harland et al. [1989] timescale, 
whereas Müller et al. [1993] used that of Kent and Grad- 
stein [1986]. We have compared the differences in the ages 
of each particular chron or anomaly, and of geological stage 
boundaries, when the two different timescales are used. 
Between 120 Ma and the Present, the corresponding mean 
age differences are generally less than 1 Myr, with a 
maximum of 1.4 Myr between 55 and 170 Ma. Between 
120 and 130 Ма (1.е., from the Valanginian to the early 
Aptian), the difference reaches 4 to 7.5 Ma, but this is a time 


Mean Age SA CN) Sọ (CE) Colatitude NC 95%CI Site Plate Rock Type Author 
3.3 —26.1 —5.1 121.7 83 4.2 519 AFR 5 Tauxe et al. [1983] 
15 20.1 61.5 82 49 4.1 222 IND 5 Pierce [1976] 
18 —26.1 —10.3 123.7 98 3.1 521 AFR 5 Tauxe et al. [1983] 
26 —26.1 —5.1 121.2 55 3l 522 AFR s Tauxe et al. [1983] 
27.5 —28 —36.3 104.7 42 44 17 AMS 5 Pierce [1976] 
28 —26.1 —5.1 122.9 100 1.5 522 АЕК 5 Tauxe et al. [1983] 
32.5 —28.5 —2.2 123.8 65 4.4 523 AFR 5 Tauxe et al. [1983] 
40.3 —28.5 —2.2 132.3 82 2.7 523 AFR 5 Tauxe et al. [1983] 
45.5 —28.5 —2.2 133.5 102 3.8 523 AFR 5 Tauxe et al. [1983] 
47 —42 73.4 102.3 52 1.8 713 IND у Vandamme and Courtillot [1990] 
54.7 —24.9 87.4 141.9 49 3 253 IND s Pierce [1976] 
57.5 —19.2 99.3 119.7 38 4.8 212 AUS s Pierce [1976] 
63.9 —29.5 3.5 135.9 257 2.6 524 AFR s Tauxe et al. [1983] 
66 =75 59 115.2 56 4.2 707 AFR у Vandamme and Courtillot [1990] 
104.5 —16.1 116.3 116.7 61 3.9 260 AUS 5 Pierce [1976] 
114.5 —23.3 111 135.9 44 4.7 263 AUS 5 Pierce [1976] 
114.5 —23.3 111 132 106 2.9 263 AUS 5 Pierce [1976] 


“SXCN), site latitude; Só(^E) site longitude; colatitude, paleomagnetic colatitude; NC, number of samples; 95%CI, confidence interval; site, number, 
plate to which datum belongs; rock type, s, sediment ; v, volcanics; author, reference. 
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Table 3. Finite Rotation Poles for Gondwana* 
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NW Africa onto South Africa 


Arabia onto South Africa 


An Аве, Ma  X(CN) (CE) w, deg Source Age, Ма ACN) o(°E) w, deg Source 
A34 84.0 0.0 0.0 0.0 4.7 32.75 22.64 —1.89 (2) 
100.0 —17.3 —172.8 5 (1) 13 32.15 22.59 —5.36 (2) 
110.0 —15.7 —174.0 .9 (1) 30 32.11 22.57 —7.36 (2) 
MO 118.7 —18.1 —173.6 L4 (1) 
MIO 130.0 29.8 15.2 1.6 (1) 
East Antarctica to South Africa India to South Africa Australia to South Africa 
An Age Ma NCN) (CE) w, deg Source ACN) (CE) w, deg Source ACN) o(°E) w, deg Source 
5 10.4 7.9 —49.5 1.5 (1) —23.7 —146.7 4.6 (1) —11.9 —130.4 6.6 (1) 
6 20.5 10.8 —48.0 2.8 (1) —29.6 —156.2 7.6 (1) —]4.5 —133.5 11.76 (1) 
13 35.5 14.1 —50.7 S5 (1) —21.4 —139.3 15.3 (1) —13.3 —130.6 20.16 (1) 
18 42.7 15.5 —44.8 7.2 (1) —23.6 —138.1 17.9 (1) —15.9 —130.8 22.76 (1) 
21 50.3 10.3 —42.9 8.8 (1) —23.7 —144.1 22.1 (1) —16.0 —127.6 23.76 (1) 
25 58.6 5.3 —39.1 10.3 (1) —23.3 —152.7 28.6 (1) —15.9 —122.7 24.16 (1) 
31 68.5 1.1 —41.6 11.8 (1) —20.5 —160.8 42.4 (1) —]11.4 —118.2 25.66 (1) 
33y 73.6 —1.9 —41.4 13.5 (1) —222 —162.3 45.7 (1) —13.5 —115.2 26.36 (1) 
330 80.2 —4.6 —39.7 16.0 (1) —23.9 —163.2 48.8 (1) —16.3 —110.5 27.16 (1) 
34 84.0 —2.6 —38.1 17,9 (1) —22.9 —161.2 51.0 (1) —14.7 — 106.3 27.86 (1) 
90.0 —3.9 —332 22.1 (1) —22.5 —156.4 53.0 (1) —15.7 —97.8 29.06 (1) 
100.0 —5.6 —29.0 29.3 (1) —21.4 —149.3 56.8 (1) —18.0 —85.0 31.86 (1) 
110.0 —6.5 —26.1 36.6 (1) —21.4 —149.3 56.8 (1) —20.5 —73.8 35.36 (1) 
MO 118.7 —7.1 —24.4 43.1 (1) —21.7 —148.4 57.0 (1) —22.0 —65.7 39.26 (1) 
MIO 130.0 —11.9 —21.0 50.4 (1) —22.3 — 143.1 61.0 (1) —28.2 — 56.9 43.56 (1) 
M16 140.0 —7.0 —26.9 50.70 (3) —21.3 —139.0 64.7 (4) —21.7 —63.8 46.9 (5) 
M21 150.0 —4.7 —29.0 52.84 (3) —21.6 —135.2 66.6 (4) —192 —63.7 49.7 (5) 
FIT 165.0 —7.8 —31.4 58.0 (5) —26.9 —133.4 67.8 (4) —22.6 —62.7 552 (5) 


Madagascar to South Africa 


An Аве, Ma ACN) (CE) w, deg Source 
MO 118.7 A —76.3 0.9 (1) 
M10 130.0 3.4 —79.2 6.9 (1) 
M16 141.9 5.40 —76.20 8.32 (6) 
M21 149.9 4.00 —71.40 11.32 (6) 
FIT 165.0 —3.41 —81.70 19.73 (7) 


“Recomputed following (1) Müller et al. [1993]; (2) Le Pichon and Gaulier [1988]; (3) Royer et al. [1988]; (4) India to South Africa motion is India to East 
Antarctica [Royer and Sandwell, 1989] plus East Antarctica to South Africa [Royer et al., 1988]; (5) Australia to South Africa motion is Australia to East 
Antarctica [Royer and Sandwell, 1989] plus East Antarctica to South Africa [Royer et al., 1988]; (6) Royer et al. [1988]; and (7) Lawyer and Scotese [1987]. 


window with no sedimentary data. From 135 Ma back to the 
earliest seafloor spreading anomalies at 160 Ma, the differ- 
ences are between 0 and 2.8 Myr, averaging 1 Myr in 
absolute value. Prior to 160 Ma, only geological age assign- 
ments of formations from which paleomagnetic data are 
available are used and not plate kinematic parameters. For 
our purpose, these differences are smaller than the uncer- 
tainties in age determinations in most paleomagnetic stud- 
ies, and often less than the error due to reassigning an age to 
an anomaly when the exact time feature used in pointing the 
anomaly is not accurately known. We have therefore used 
ages as given by the original authors, regardless of the 
timescale they used. However, as mentioned above in the 
discussion of the new paleomagnetic database, we have 
scaled the ages of all DSDP/ODP sites using the Harland et 
al. [1989] timescale. 


4. New Synthesized APWPs 
4.1. African APWP: Old and New 

[24] All data and sites are shown in Figures la and 1b in 
(South) African coordinates, together with average poles 
(Table 4) calculated every 10 Myr with a 20 Myr window, 
using a method derived from McFadden and McElhinny 
[1988] in order to combine poles and small circles (this is 


done by iteratively finding a direction such that the Fisher 
average of pole directions and the points on the small circle 
which are closest to this direction are statistically the same). 

[25] In order to give some feeling for the quality of 
individual means, four examples showing the poles and 
sites on an appropriate plate reconstruction are given for 20, 
60, 90, and 200 Myr in Figure 2. Thirty-nine individual data 
are available for the 20 Myr reconstruction (Figure 2a), with 
good longitudinal and latitudinal coverage. The correspond- 
ing pole distribution is Fisherian and the 95% confidence 
interval is less than 3°. Data come from most major plates. 
The 60 Myr reconstruction (Figure 2b) involves 46 data; the 
distribution is approximately Fisherian with a 95% confi- 
dence interval of 3°, although several poles streak E-W, east 
of (65°N, 270°E). There are no data from South America. 
The 90 Myr reconstruction (Figure 3a) contains only 13 
data with a 95% confidence interval of 2.9°, but geo- 
graphical coverage is quite good, with all major landmasses 
except Antarctica providing compatible data. There are still 
20 data available for the 200 Myr reconstruction (Figure 3b), 
again with most plates (except South America) in their 
Pangea configuration represented, a good Fisherian pole 
distribution and a 95% confidence interval of about 4°. 
Altogether, the Ao; values range from 2 to 7°, with only a 
slight degradation as one goes back in time (Figure 1d). The 
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Table 4. Master Apparent Polar Wander Path for the Past 200 Myr Calculated for a 20 Myr Sliding Window Every 10 Myr? 


South Africa 
Window Age М X(N) ф (Е) As К X(N) 6 (E) 


South America India 
Ags K 


Australia 
ACN) 6 CE) 4o5 K A(N) ФСЕ) Aos K 


0 3.1 30 865 1716 26 105.3 86.3 172 
10 8.3 54 860 160.8 20 958 854 162.5 
20 189 38 854 151.9 2.7 75.9 840 154.8 
30 295 23 85.1 1622 38 656 82.8 158.1 
40 40.0 24 843 1724 3.3 82.1 81.33 162.4 
50 52.2 31 847 1747 34 601 80.9 164.4 
60 597 45 847 2176 28 574 81.1 190.5 
70 67.3 34 83.8 2416 3.2 603 803 204.3 
80 77.9 14 847 2758 60 451 814 206.1 
90 90.0 13 85.0 3208 53 617 822 202.1 
100 97.6 12 85.5 9.7 68 42.0 817 180.1 
110 113.6 17 777 20.2 4.1 76.7 80.0 183.6 
120 119.1 20 764 173 2.3 209.6 78.2 189.4 
130 126.4 14 753 140 32 1545 75.8 192.9 
140 136.8 7 724 5.8 65 87.4 73.8 197.6 
150 151.6 10 67.0 26.6 6.8 50.8 750 159.9 
160 1623 15 629 31.6 5.0 585 72.5 144.0 
170 173.4 21 56.9 393 6.0 28.8 697 112.5 
180 178.8 18 53.2 45.7 54 413 65.5 95.9 
190 189.7 23 545 45.0 4.2 529 65.3 98.4 
200 196.7 19 58.2 46.9 43 61.6 63.2 106.0 


2.6 105.1 86.1 1748 26 1052 861 1748 2.6 105.2 
20 944 85.0 168.1 2.0 94.2 85.0 1681 2.0 94.2 
2.7 76.2 833 1642 2.7 75.6 833 1642 2.7 75.6 
3.8 662 81.5 1692 3.8 654 81.5 1692 3.8 65.4 
3:3 81.9 79.5 1744 3.2 85.5 791 17531 32 87.5 
34 594 77.9 1793 3.4 58.2 76.3 1780 3.4 59.7 
2.99 561 759 196.8 2.9 547 744 1913 2.8 57.2 
32 61.0 742 2048 3.2 61.2 72.9 197.9 32 61.7 
5.9 47.2 747 2074 5.9 47.0 73.1 2022 60 45.7 
5.2 652 75.5 2074 5.1 65.9 729 203.8 51 67.2 
6.7 43.0 76.6 195.8 6.7 43.1 73.0 1944 6.7 43.0 
42 748 751 193.8 4.2 75.3 713 1945 42 75.0 
24 182.9 73.1 193.99 24 1843 693 196.0 24 183.7 
28 205.5 70.6 193.0 2.8 205.6 66.8 196.5 2.8 205.8 
6.0 103.2 68.3 1942 60 1034 646 1984 60 103.4 
6.66 543 73.6 1677 6.6 54.5 687 1752 66 544 
5.0 597 737 1497 5.0 597 684 1613 50 597 
6.7 236 755 1101 6.7 23.6 709 131.8 67 23.6 
5.66 39.7 73.0 834 5.6 397 69.9 1092 5.6 39.7 
4.2 52.9 726 868 42 529 693 1118 42 52.9 
4.3 61.6 69.8 956 43 61.6 66.1 1177 43 61.6 


Antarctica Europe North America Greenland 
Window Age N ACN) ф (Е) As К ACN) (CE) As K ACN ФСБ) As К (№) Ф (Е) As K 
0 3.1 30 865 1716 26 105.3 86.3 172 2.6 105.1 86.1 1748 26 1052 861 1748 2.6 105.2 
10 8.3 54 860 1608 2.0 958 854 1625 20 944 85.0 1681 2.0 942 85.0 1681 20 94.2 


20 18.9 38 854 1519 27 759 840 154.8 
30 295 23 85.1 1622 38 656 828 158.1 
40 40.0 24 843 1724 3.3 82.1 81.3 162.4 
50 52.2 31 847 1747 34 601 80.9 164.4 
60 59.7 45 847 2176 28 574 81.1 190.5 
70 67.3 34 83.8 2416 32 603 80.3 204.3 
80 77.9 14 $847 2758 60 451 814 206.1 
90 90.0 13 85.0 3208 53 617 822 202.1 
100 97.6 (12 85.5 9.7 68 420 817 180.1 
110 1136 17 777 20.2 41 76.7 800 183.6 
120 1191 20 764 17.3 23 2096 78.2 189.4 
130 1264 14 75.3 14.0 32 154.5 75.8 192.9 
140 136.8 7 72.4 5.8 65 874 73.8 197.6 
150 1516 10 67.0 26.6 68 50.8 75.0 159.9 
160 162.3 15 629 31.6 5.0 58.5 72.5 144.0 
170 173.4 21 569 39.3 6.0 28.8 697 112.5 
180 1788 18 532 457 54 41.3 65.5 95.9 
190 189.7 23 54.5 45.0 4.2 52.9 65.3 98.4 
200 196.7 19 58.2 46.9 43 61.66 63.2 106.0 


2.7] 76.2 83.3 1642 2.7 75.6 833 1642 27 75.6 
3.8 662 81.5 1692 3.8 654 81.5 1692 3.8 654 
3.3 81.9 79.5 1744 32 855 79.1 1751 32 87.5 
34 594 77.9 1793 34 582 76.3 1780 3.4 59.7 
2.99 56.1 75.9 1968 2.9 547 744 1913 2.8 57.2 
32 610 742 2048 3.2 61.2 729 1979 3.2 61.7 
59 472 747 2074 5.9 47.0 731 2022 60 45.7 
5.2 652 75.5 2074 5.1 65.9 729 203.8 5.1 67.2 
6.7 43.0 76.6 195.8 6.7 43.1 73.00 1944 67 43.0 
42 748 75.1 193.8 4.2 75.3 713 1945 42 75.0 
24 182.9 73.1 193.9 24 1843 693 196.0 24 1837 
28 205.5 70.6 193.0 2.8 205.6 66.8 196.5 2.8 205.8 
6.0 103.2 68.3 194.2 60 1034 646 1984 60 103.4 
6.66 543 73.6 1677 6.6 54.5 687 1752 66 544 
5.0 597 737 1497 5.0 59.7 684 161.3 5.0 59.7 
6.7 23.6 75.5 1101 6.7 23.6 70.9 131.8 67 23.6 
5.6 39.7 73.0 834 5.6 39.7 69.9 1092 5.6 39.7 
42 52.9 726 868 4.2 529 693 1118 42 52.9 
43 61.6 69.8 956 43 61.6 66.1 1177 43 61.6 


“Window, age of the center of window; age, mean age computed from the data; N, number of studies; ^, ф, latitude and longitude of mean VGP; Ао» 


uncertainty at the 9596 confidence level; K, Fisher's precision parameter. 


overall mean Ags is 2.9 + 0.8? (16) for the new synthesized 
APWP, significantly less than the value of 4.1 + 1.2? found 
by BC91; the decrease (4.1/2.9 = 1.4) is not much less than 
the maximum value (4/2.2 = 1.5) expected from increasing 
the number of data. 

[26] There are enough data that a significant synthetic 
APWP can be calculated with a time resolution increased 
by a factor of 2, i.e., one average every 5 Myr with a 10 
Myr time window (Table 5). This results in a mean Ags of 
4.2 + 1.9°, i.e., similar to BC91 with a twofold increase in 
time resolution. 

[27] We have plotted in Figure 4 the angular test param- 
eter үсү, [from McFadden and McElhinny, 1990], deter- 
mined for pairs of poles consisting of the calculated mean 
pole from each plate versus the corresponding overall mean 
(synthesized) pole for the same 20 Myr time window, as a 
function of age. Pole pairs are not significantly different at 


the 95% probability level when the test parameter is 
positive. Four apparently discrepant poles deserve some 
comment. Two poles appear to be significantly different 
from the coeval synthesized pole: those for NAM in the 0— 
10 Myr window and EUR at 70 Ma (60—80 Myr window). 
For the NAM pole, there are only 6 data and the discrep- 
ancy could for instance be due to a slightly larger quadru- 
pole term, on the order of 10% of the axial dipole (for the 
last few million years, ranges of estimates are between 3 and 
8% [Schneider and Kent, 1988; McElhinny et al., 1996; 
Johnson and Constable, 1997; Carlut and Courtillot, 
1998]). The pole from EUR is based on a combination of 
only three studies, all being to some extent questionable. 
One study comes from the Aix-en-Provence series (south of 
France), in which rotations have subsequently been 
revealed. The two other studies are from the Antrim and 
Mull lavas of the British Tertiary Igneous Province. Their 
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b (60.0 Ma) 


Figure 2. Plate reconstructions at (a) 20 and (b) 60 Ma in rectilinear projection, together with 
transferred poles in equal-area projection (the South African plate is kept fixed in longitude). The 
transferred poles and site locations are also shown as open circles and solid stars, respectively, in 
the projections. Insets: polar view of the synthetic African APWP (small linked dots) with all data for the 


appropriate time window (larger open circles). 


ages in the database (around 70 Ma) are now known to be 
too old by some 10 Myr when compared with more recent 
estimates based on a combination of new radiometric ages 
and the magnetic polarity timescale [Saunders et al., 1997]. 
Indeed, a large part of these lavas erupted during reversed 
chron C26r dated at 59 + 1 Myr. Such an error in the 


assigned age may easily account for the observed discrep- 
ancy. For two poles, there is a problem of rounding the age 
at a time of significant change. Such is the case for the 50— 
60 Ma Australian pole and for the 160—170 Ma Greenland 
pole. The effect of the slight differences in ages of observed 
versus predicted poles illustrates difficulties linked to the 
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b (200.0 Ma) 


Figure 3. Plate reconstructions at (a) 90 and (b) 200 Ma. Same legend as Figure 2. 


sparseness of the database and can lead to misleading 
tectonic or field-geometry interpretations. For instance the 
calculated 54.4 + 5.2 Myr mean pole from AUS data only 
can be compared either with either the 55 or 59 Myr overall 
mean (synthesized) pole, in one case with a discrepancy, in 
the other with none. A similar case is also shown for GRE at 
160—170 Ma. Note again (as in Figure 1d) that there is no 
sign of degradation with age. 

[28] In Figure 5, we compare three synthesized APWPs in 
(south) African coordinates: triangles are for the original 
BC91 curve, dots for the new ВСО1 curve, and stars for an 


APWP computed with the BC91 database but with the 
updated kinematics used in the present paper (ВС91”). The 
0—140 Ma and 120—200 Ma segments are shown separately 
in Figure 5 to avoid too much overlap and lack of legibility. We 
first see that improvement in the kinematic models results in 
insignificant changes from the Present back to 160 Ma. On the 
other hand, the 170 (actually 175) and 180 (actually 178—179) 
Ma poles are significantly displaced. This is essentially due to 
a different Gondwana fit (see section 4, APWP for Africa). 
[29] Now, comparing BC91 and BCOI, we see that the 
increased database has resulted in a systematic ~2° shift of 
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Table 5. Synthetic Apparent Polar Wander Path for the Past 200 Myr Calculated for a 10 Myr Sliding Window Every 5 Муг“ 


South Africa South America India Australia 
Window Age N (CN) ọ (Е) Aos K ACN) CE) Aos K A(N) ọ CE) 45 K ACN) ФСЕ) 45 K 
0 2. 25 866 1822 3.0 962 868 1763 3.0 96.0 87.0 1971 3.0 957 875 2040 3.0 951 
5 3.1 30 862 176.9 2.6 105.5 86.5 168.0 2.6 1048 86.8 197.1 2.6 106.2 87.5 2076 2.6 105.1 
10 11.9 21 85.0 170.7 3.1 107.4 85.5 139.9 3.1 107.3 85.8 231.1 3.1 104.0 85.9 277.2 3.2 100.9 
15 14.8 24 842 1709 3.2 85.1 84.8 1367 32 842 849 231.2 32 868 849 2793 32 89.0 
20 196 16 817 165.7 45 693 82.1 1318 46 664 834 2219 45 70.0 841 273.0 45 69.3 
25 260 14 834 1845 54 558 844 1261 53 592 81.3 2647 57 508 786 294.7 57 50.8 
30 302 13 79.0 2077 54 63.3 £835 163.8 53 64.7 749 263.6 50 723 73.0 2828 50 72.4 
35 33.7 12 78.9 2017 45 96.1 829 1521 47 894 75.1 2663 45 966 72.9 2853 45 96.5 
40 39.0 8 773 1916 72 63.3 800 1395 7.3 620 749 2644 7.2 637 732 2837 71 64.7 
45 464 12 791 2006 52 730 820 1265 53 726 706 277.8 55 67.1 702 2946 54 69.1 
50 499 17 78.7 2041 42 745 822 1277 42 755 677 2785 45 66.1 69.5 295.7 4.1 76.4 
55 550 22 767 2147 42 56.1 827 1453 41 57.5 607 277.7 48 43.6 672 2942 41 58.0 
60 60.7 24 73.9 2240 43 48.7 824 1685 43 48.9 508 2775 46 425 643 292.7 43 48.8 
65 64.0 24 71.6 2342 36 67.8 820 1937 3.6 68.3 42.4 2794 4.0 56.1 60.6 293.1 3.6 67.9 
70 683 15 711 2374 47 672 823 1990 48 65.0 368 2810 45 745 595 2954 48 65.8 
75 754 10 72.5 2324 74 449 836 1706 7.2 464 33.6 2821 7.7 41.2 621 3006 73 45.4 
80 791 9 683 251.9 70 56.1 823 2284 69 57.5 257 287.5 71 542 548 3035 69 58.1 
85 834 5 669 2531 98 617 82.5 2295 102 57.7 226 2891 94 672 544 3054 10.0 59.4 
90 914 8 66.7 2427 34 2598 848 175.8 3.6 2326 21.9 289.2 3.9 2000 59.9 309.4 3.9 203.5 
95 941 8 665 2460 66 722 862 1785 65 740 20.5 2917 7.0 63.9 594 3136 66 711 
100 100.2 6 688 2541 11.0 396 87.8 292 11.5 35.9 19.9 298.0 10.9 40.1 592 327.5 11.9 33.8 
105 1046 4 62.9 251.5 256 154 880 179.8 246 16.6 149 295.0 256 154 586 322.7 243 17.0 
110 110.8 7 58.0 2617 83 57.1 85.1 2642 8.7 52.2 8.3 297.5 83 56.9 53.8 3262 8.9 49.9 
115 1164 13 563 2609 2.7 2355 864 2612 24 3139 68 2972 2.7 2405 541 330.5 23 342.0 
120 120.1 13 51.8 260.9 3.2 1639 83.9 2385 31 1782 2.4 296.1 32 166.2 52.9 327.9 2.9 204.6 
125 122.8 10 506 2607 34 1975 834 2330 34 1977 0.9 296.4 3.7 1739 53.1 3286 3.8 166.4 
130 130.9 4 495 265.8 74 1569 82.6 2526 72 1630 —3 301.5 8.0 1323 50.2 3356 382 127.3 
135 1343 5 495 2646 56 1892 83.2 2465 5.5 1926 —35 3016 59 1675 49.6 3354 62 152.3 
140 1393 5 436 2659 81 90.1 77.5 2382 81 90.1 —10 3007 79 93.8 1441 3285 82 884 
145 1420 3 403 263.2 11.9 1081 747 2272 11.9 1081 -—12.4 297.9 117 1117 43.1 3232 12.0 107.0 
150 153.9 4 591 2557 63 2162 85.8 616 63 2162 5.7 3083 6.0 2319 506 3524 61 231.1 
155 155.7 7 58.2 258.9 3.8 253.8 87.3 396 3.8 253.8 4.1 3098 3.9 2446 485 351.8 3.9 245.7 
160 1589 5 53.8 257.9 71 1176 88.2 1896 7.1 117.6 1.1 3079 64 145.8 479 3466 64 144.1 
165 168.0 8 521 2605 98 33.7 866 2285 98 337 —07 3107 98 337 448 3473 98 33.7 
170 170.2 13 565 259.1 68 378 88.9 41.9 68 37.8 3.4 3124 68 378 462 353.3 68 37.8 
175 176.7 13 65.7 260.8 72 33.8 79.7 322 7.2 33.8 10.6 3183 7.2 33.8 45.5 6.6 7.2 33.8 
180 178.4 12 668 2639 76 334 784 266 7.6 33.4 10.9 3199 76 334 442 8.0 7.6 33.4 
185 1824 9 695 2614 76 465 75.9 326 7.6 46.5 13.6 3206 76 465 449 1.9 7.6 46.5 
190 191.3 13 629 2589 55 584 82.5 37.0 5.5 58.4 8.8 3160 5.5 584 464 2.6 5.5 58.4 
195 1944 14 617 2579 49 679 8327 416 49 67.9 8.0 3149 49 679 46.9 0.8 49 67.9 
200 1989 8 63.9 2446 5.7 949 788 70.06 57 94.9 13.7 3117 57 949 52.6 5.8 57 949 
Antarctica Europe North America Greenland 
Window Age N CN) CE) Aos К ACN Фф (Б) dos К ACN) CE) Aos К A(N _ CE) As К 

0 2. 25 868 1788 3.0 960 867 1787 3.0 96.1 865 1807 3.0 962 865 1807 3.0 96.2 
5 31 30 865 171.6 2.6 105.3 863 1720 26 105. 86.1 174.8 2.6 1052 86.1 174.8 2.6 1052 
10 11.9 21 85.9 151.8 31 1073 85.0 1557 3.1 1076 846 1644 3.1 1077 846 1644 31 107.7 
15 148 24 85.2 151.0 32 85.1 84.2 1549 32 843 836 163.0 32 842 836 1630 32 842 
20 196 16 828 1467 45 679 814 149.7 45 67.8 810 1562 45 68.3 81.0 156.2 45 68.3 
25 260 14 858 1521 54 575 83.8 1532 5.3 58.7 828 1657 53 579 82.8 1657 53 57.9 
30 302 13 83.5 197.0 52 66.1 816 1834 53 644 796 1879 54 63.2 79.6 187.9 54 63.2 
35 33.7 12 83.6 1853 46 916 812 1734 46 904 79.3 1804 46 924 793 1804 46 92.4 
40 39.0 8 81.8 1672 73 62.7 788 1602 7.3 622 773 1677 73 624 772 1680 73 61.8 
45 464 12 848 1549 53 71.2 81.1 1504 5.3 727 79.6 167.9 52 735 78.7 169.9 51 76.3 
50 499 17 85.1 1564 42 76.0 813 151.9 42 75.3 79.3 1703 42 744 780 171.9 42 73.9 
55 550 22 853 181.7 41 58.1 814 1683 42 570 77.9 1834 42 55.9 76.1 180.9 4.2 56.9 
60 60.7 24 843 2115 43 489 805 1889 43 488 754 1955 43 48.7 73.8 1893 43 48.9 
65 64.0 24 827 2406 3.6 68.1 79.8 209.5 3.6 68.0 73.5 207.3 3.6 67.9 72.6 2001 3.6 684 
70 683 15 82.8 251.7 48 652 800 2132 48 645 734 209.7 48 649 72.5 2023 48 65.8 
75 754 10 867 245.7 73 45.8 81.3 188.6 72 469 757 1976 72 468 737 192.0 73 46.0 
80 79.1 9 80.9 292.0 69 58.1 81.0 2325 69 57.0 73.5 2214 69 570 725 2160 69 568 
85 834 5 81.0 298.2 102 569 81.1 2305 103 563 73.5 220.5 102 567 721 216.1 102 57.6 
90 914 8 876 3249 3.9 2034 80.8 185.8 3.5 2456 75.3 1964 3.5 2471 71.9 195.5 3.5 2517 
95 941 8 863 3564 66 722 818 1834 64 765 764 197.7 64 763 73.00 195.8 64 76.3 
100 1002 6 80.3 32.6 12.0 334 85.1 144.5 11.5 36.3 817 198.3 11.5 363 78.1 1904 11.4 36.5 
105 1046 4 824 23.7 242 17.1 812 173.5 247 165 768 191.5 246 165 73.0 1913 247 16.5 
110 110.8 7 7841 9.2 89 49] 821 2038 87 519 754 208.8 87 517 72.5 2061 87 51.8 
115 1164 13 76.5 19.6 22 23465 79.7 1862 2.5 2784 747 1946 25 2813 70.9 1955 2.5 2803 
120 120.1 13 77.0 11.9 2.9 2059 765 193.5 3.1 1785 71.2 1942 31 1786 674 1972 31 1783 
125 122.8 10 76.8 142 3.8 1657 75.7 1924 3.5 193.5 705 1925 3.5 193.6 667 1961 3.5 193.8 
130 1309 4 719 18.8 7.8 1397 774 1952 7.0 1728 71.9 1965 7.0 173.8 68.2 1987 7.0 1743 
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Antarctica Europe North America Greenland 
Window Age N A(N) ó(E) dos K  X(N) ФОБ) 45 К ACN) (CE) 4o К ACN) ФСЕ As K 
135 1343 5 718 185 60 165.66 76.8 190. 5.5 1945 71.8 1923 55 195.1 67.9 1953 5.5 195.7 
140 1393 5 70.9 358 80 923 726 2031 81 900 66.7 1972 81 90.0 63.2 2016 81 90.0 
145 142 3 720 3468 119 109.1 68.7 201.9 11.7 113 63.1 193.0 11.7 1112 59.4 199.8 11.7 111.2 
150 153.9 4 63.9 41.5 61 2310 732 1343 61 2272 757 141.6 61 2269 704 1552 6.1 226.6 
155 155.7 7 62.9 37.2 3.9 245.6 743 1374 3.5 2962 76.2 1468 3.5 2955 709 1591 3.5 295.2 
160 158.9 5 65.2 31.1 64 1441 719 1497 71 1179 724 1541 71 1177 672 1648 71 1178 
165 168 8 62.8 265 98 337 706 1495 9.7 340 71.2 1519 97 340 659 1632 97 34.0 
170 170.2 13 60.6 348 68 378 70.1 1299 7.7 304 735 131.5 7.7 304 684 1474 77 30.4 
175 176.7 13 53.2 45.1 72 33.8 660 964 74 324 734 845 74 324 703 1104 74 324 
180 178.4 12 51.6 446 76 334 660 917 76 334 73.8 77.9 76 334 711 105.2 7.6 33.4 
185 182.4 9 50.1 485 76 465 637 877 76 465 71.7 71.2 76 465 69.6 97.9 7.6 46.5 
190 191.3 13 55.9 43.1 55 584 659 1025 55 584 72.8 928 55 584 692 1168 5.5 584 
195 1944 14 57.2 42.3 49 67.9 65.8 1058 49 67.9 72.5 972 49 679 68.6 1202 4.9 67.9 
200 198.9 8 580 544 57 949 594 1035 5.7 949 664 90.5 5.7 949 630 1119 57 949 


* Same as Table 4. 


the 10 to 60 Ma poles, and ~4° shift of the 80 to 110 Ma 
poles. The 120 and 130 Ma poles agree, but the two paths 
diverge prior to 140 Ma, with a maximum difference of 14° 
at 170 Ma. Actually, when one takes into account the Aos 
uncertainties, only the 173—175 and 178 Ma poles fail to 
intersect. Altogether, the new ВСО1 African APWP returns 
on itself rather than forming an open loop in the 90—200 Ma 
interval. 


4.2. Checking the Geocentric Axial Dipole Hypothesis 


[30] The comparison between the 10 Myr and 5 Myr 
resolution APWPs is best described in the case of the plate 
with the fastest polar wander rates, i.e., India. This is done 
in a later section. However, before comparing the new and 
old versions of the synthetic APWP for other plates, we 
address the question of dipolarity of the past geomagnetic 
field. As recalled above, a significant quadrupole, on the 
order of a few percent of the axial dipole, has been 
identified in the mean geomagnetic field for the last 5 
Myr [e.g., Johnson and Constable, 1997; Carlut and 
Courtillot, 1998]. In Figure 6, we have plotted all the poles 
from our database referred to a common site longitude. In 
each 20 Myr time window, all sites are related to a common 
longitude (taken to be 0), and all such frames are stacked 
from 0 to 200 Ma. When a mean pole is computed in each 
time window, a jagged path results (Figure 7). Poles tend to 
follow an erratic course, yet remaining most of the time in 
the hemisphere opposite to the one centered on the reference 
meridian. Of course, the angular distance from the pole is 
small (on the order of 2° on average) and the uncertainties 
are such that most 95% confidence intervals include the 
geographical pole. When a grand 200 Myr average is 
computed, the mean pole lies at X = 88.6°N, ф = 176.5°E 
(Ags = 1.2°), i.e., on the far side of the reference meridian, 
1.4 + 1.2° away from the geographical pole. If taken at face 
value, this would imply a far-sided effect due to a persistent 
quadrupole on the order of 3% (42%) of the dipole, i.e., 
somewhat smaller than the 5 Myr average. This translates to 
a maximum effect of 3° on inclination at the equator, and 
translates into a paleolatitude error of 1.5° at most when the 
axial dipole hypothesis is used. We conclude that there are 
indications in favor of a persistent quadrupole on the 200 
Myr timescale, with an amplitude on the order of half the 
one obtained for the last 5 Myr. However, these values are 


not distinguishable from zero in any individual 10 or 20 
Myr time frame (except at 52 and 59 Myr), and in any case 
the effect is small and on the order of contributions from 
other sources of paleomagnetic uncertainty. With the reso- 
lution we are able to achieve, these quadrupolar effects 
remain negligible and the geocentric axial dipole (GAD) 
hypothesis is clearly satisfactory on the timescales of 
interest to our study (10 to 200 Myr). 


4.3. Description of Individual Plate APWPs 


4.3.1. Africa 

[31] We have already illustrated the main characteristics 
of the BCO1 synthetic APWP in South African coordinates 
(using the 20 Myr sliding window) (Figures | and 5). The 
path can be described in terms of a succession of track 
segments and standstills, i.e., of episodic polar wander [e.g., 
Briden, 1967; Gordon and Cox, 1980; Cox and Hart, 1986; 
Irving and Irving, 1982]. Poles tend to cluster at 10—20 Ma, 
90—100 Ma, 180—190 Ma and a sharp hairpin occurs at 
140 Ma. These relative standstills and cusps separate tracks 
with rather fast and regular polar wander: 20—90 Ma, 100— 
140 Ma, 140—180 Ma. 

4.3.2. Europe and North America 

[32] The synthetic APWPs in European and North Amer- 
ican coordinates share the same major features, related 
through the slow opening of the North Atlantic Ocean at 
the end of the Cretaceous (Figure 8a). A track from 10 to 50 
Ma is preceded by a more complex path between 130—140 
and 50 Ma, with a general change in trend. A previously 
undetected loop occurs between 50 and 110 Ma. It can in a 
way be taken to represent a standstill, but some of the 
features in the loop appear to be resolved by individual 
mean poles. In particular, the 52 and 59 Ma averages are 
statistically distinct. A problem comes from the fact that 
much of the loop overlaps in time with the Cretaceous Long 
Normal Superchron, during which oceanic spreading rates 
are assumed to be constant for lack of observable chrons. 
This introduces a decrease in the temporal resolution of the 
kinematic models which are used to transfer poles. 

[33] There are distinct cusps at 140 and 180—190 Ma, as 
was the case for the African path. The overall shape of the 
Eurasian APWP between, say, 50 and 150 Ma confirms the 
loop which was one of the new, previously ill-recognized 
features uncovered in BC91. The loop appears a bit smaller 
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Figure 4. The McFadden and McElhinny [1990] test applied to each calculated mean pole for each 
plate versus the corresponding overall mean (synthetic) pole for the same 20 Myr time window as a 
function of age (the vertical axis corresponds to the calculated test parameter у. ү). If negative (gray 
shaded zone), pole pairs are considered distinct at the 9596 probability level. Two poles (AUS near 55 Ma 
and GRE near 165 Ma) fall in between rounded ages and can be assigned to either one of two synthetic 
poles (at 50 and 60 Ma and 160 and 170 Ma, respectively); hence two test points linked by a straight line 


segment are shown. 


and may have a more complex structure than previously 
recognized. The timing of the so-called Cretaceous standstill 
[Westphal et al., 1986; Besse and Courtillot, 1991] appears to 
be somewhat later: 60—110 Ma rather than 70—130 Ma. This 
loop was a significant feature used to analyze paleomagnetic 
data from mobile zones and large continental blocks in Asia 
[e.g., Besse and Courtillot, 1988, 1991; Enkin et al., 1992]. 
Consequences in the changes in shape ofthe Eurasian APWP 
for such analyses are not explored further in this paper. 

[34] Because most data on which the Jurassic segment of 
the NAM APWP could be based carried some amount of 
uncertainty (tectonic rotations for data from the SW United 
States, remagnetization on intrusions from the NE United 
States), Courtillot et al. [1994] proposed to transfer their 
available data from other continents to generate synthetic 
NAM APWPs. They discussed the effects of data selection 
and kinematic reconstructions. They concluded that the high 
latitude APWP of /rving and Irving [1982] was to a large 
extent vindicated, though with much reduced uncertainties at 
the 95% confidence level. On the other hand, the lower 
latitude paleomagnetic Euler pole-based path of May and 
Butler [1986], though on the edge ofthe confidence intervals, 
was not supported by the synthesized path. It is interesting to 
compare the path for North America from the present study 
with the earlier attempts by Courtillot et al. [1994]. Basically, 
the new path (Figure 8b and Table 4) vindicates the con- 
clusions based on BC91. The 130—200 Ma segment of the 
NAM APWP remains at rather constant latitudes on the order 
of 70—75?N, with the mean Ags reduced from 7° to 5°. 
Reliable poles are now obtained at 140 and 170 Ma, due to 
the increased number of data, which were not sufficient in 
BC9]. The main difference lies with the 200 Ma pole. Van 


der Voo [1990] emphasized that a pole transfer using the 
Bullard et al. [1965] fit led to North American and European 
APWP segments in better agreement during the Paleozoic 
and part of the Mesozoic. The problem is to understand when 
and how the change in configuration between the reconstruc- 
tion based on the oldest seafloor data [see Royer et al., 1992] 
on the one hand, and on the Bullard et al. [1965] fit on the 
other, was achieved. Our European database allows one to 
compute separate 185 and 190 Ma mean poles with three and 
five studies, respectively. The comparison between these two 
poles and a version of our master APWP from which Euro- 
pean data have been removed is inconclusive at the 95% 
level, leaving this question still open, as already concluded 
by Van der Voo [1993]. 
4.3.3. South America 

[35] As is well known, the SAM APWP is remarkable in 
that it displays very little polar wander (Figure 9a). Most 
poles have latitudes higher than 80°: South America has 
basically remained in the same position with respect to the 
geographical poles (1.е., spin axis) in the last 200 Myr. 
However, the path is ordered and shows small but resolv- 
able features, such as a track from 200 to 180 Ma, a change 
in direction near 180 Ma, a track from 180 to 130—140 Ma 
and little motion since (with a real standstill between 50 Ma 
and the Present). 
4.3.4. Australia, India, and Antarctica 

[36] It is particularly informative to draw the three 
APWPs on Figure 9b. In BC91, the Indian APWP was 
shown only back to 120 Ma, and only two paleomagnetic 
poles for that plate (Deccan and Rajmahal Traps) were 
available; Australia and Antarctica had not been included 
in our previous analysis. 
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Figure 5. Comparison of three synthesized APWPs in South African coordinates with their associated 
95% ellipses of confidence: triangles are for BC91 [Besse and Courtillot, 1991], dots for BCO1 (this 
paper), and stars for BC91" (BC91 data transferred using the kinematics used in the present paper). (a) 
From Present to 140 Ma and (b) from 120 to 200 Ma. The ellipses of confidence for the BC91' path are 
not figured for more clarity. In Figure 5b, the age numbers in italic refers to the BC91’path. Ages are in 


Ma. Equal-area projection. 


[37] The Indian path is the most elongated one, corre- 
sponding to the fastest average apparent polar wander veloc- 
ity. Three phases (tracks) are clearly recognized from 180 to 
130, 130 to 80—70 and 80—70 Ma to the Present. These are 
separated by a standstill (1.е., poles are in the joint intersection 
of their 95% confidence intervals) at 200—170 Ma, another 
one at 140—110 Ma and a directional change at 80—70 Ma. 
The classical rotation and northward motion of India are 
clearly outlined; the fast velocity (150 km/Myr) peaking 
between 70 and 50 Ma, and the slowdown after the India- 
Asia collision at 50 Ma are also seen. The 66 Ma Deccan 
Traps pole [Vandamme et al., 1991], which is based on a very 
large number of data already used in our Indian selection of 


poles (and for that reason, not used in constructing the 
synthetic path) and very accurately dated falls only 1? away 
from the 67 Ma synthetic pole (Table 4), attesting to the 
validity of the data selection process, to the quality of the 
selected paleomagnetic data and plate kinematic parameters. 

[38] Because of its fast polar wander rate, the synthetic 
APWP in Indian plate coordinates is the one for which the 
influence of improving time resolution from 20 Myr to 10 
Myr moving averages is best displayed. This is shown in 
Figure 10, where only the 95% confidence circles for the 
10 Myr APWP have been included for clarity. The mean 
poles are labeled with the mean age of the data they are 
actually based on (Tables 4 and 5). The excellent agree- 
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Figure 6. All poles from our database referred to a 
common site longitude. Open star and 95% confidence 
circle: grand 200 Ma average. Inclination alone allows only 
the determination of small circles of equal paleo-colatitude 
on which the pole must lie. 


ment between the two curves is obvious. Of course, some 
features are more readily apparent in the higher resolution 
curve, such as the slowdown at 55 Ma, compatible with an 
early collision age for India versus Asia [e.g., Patriat and 
Achache, 1984; Jaeger et al., 1989], or the sharper reorien- 
tation in the track at 75—80 Ma. The 10 Myr resolution pole 
at 100 Ma (nominal age 100.2 + 4.1 Ma, with 6 data) and the 
corresponding 20 Ma resolution pole (nominal age 97.6 + 5.8 
Ma, with 12 data) are both well within the intersection of 
their respective 95% confidence circles (10.9 and 7.3°). The 
small loop in the 10 Myr APWP based on four mean poles 
between 130 and 142 Ma is not resolved given Ags between 
6 and 12?; all are within less than 5? from the 20 Ma 
resolution pole at 140 Ma (nominal age 136.8 + 5.4 Ma) 
which has an Ags of 6.4°. At the older end of the curves, the 
10 Myr resolution poles at 194 and 199 Ma and the 20 Myr 
resolution pole at 196 Ma are fully compatible. 

[39] In conclusion, the 20 Myr resolution curve captures 
all the essential features which can be found in the 10 Myr 
curve, and it is not expected that limiting our discussion to 
the former would run the risk of having missed any 
significant feature. The high-resolution curve also confirms 
the validity of the rough classification of APW periods into 
tracks and standstills. 

[40] The Australian APWP (Figure 9b) displays almost 
the same track times and orientations, standstills and track 
changes as that for India until 70 Ma. Motion since 70 Ma is 
similar to that indicated by India (1.е., mostly northward 
motion follows after counterclockwise rotation) except it is 
three times slower between 80 and 50 Ma, and then 
accelerates a bit from 40 Ma onwards, in relation with the 
southward jump of the eastern branch of the South Indian 
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Ocean ridge. After that time, Australia and India belong to 
the same plate. 

[41] The Antarctic APWP (Figure 9b) again displays 
similar features prior to 70 Ma, with two standstills and 
two tracks, but since then, apparent polar wander is almost 
negligible, much like South America. The shapes of the 
three APWPs shown 1n Figure 9b are easily related one to 
the other as a function of the opening of ocean basins 
between plate pairs. The relative orientations of the relative 
velocity vectors and the respective rates of polar wander and 
plate motions combine in such a way that the respective 
influences of and relationships between the various param- 
eters are readily apparent. 

4.3.5. АП Paths 

[42] The succession of tracks, standstills, loops, and 
directional changes for all APWPs shows similarities and 
differences between the various continents. All paths begin 
with slow motion or standstill between 200 and 170—180 
Ma and then turn into tracks which last until 140 Ma. A 
major re-organization occurs at 140 Ma, with AFR, EUR, 
NAM and SAM changing track direction, and IND, AUS, 
ANT entering a standstill. IND, AUS, ANT and SAM 
undergo a change at 110 Ma, but the other three plates do 
not. AFR, EUR, and NAM have a track ending near 100 
Ma, and enter a standstill or loop. A 80 Myr change is 
common to IND, AUS and ANT; one at 40—50 Ma to AFR, 
EUR and NAM. Altogether, the 180—140 Ma track is the 
main feature which stands out to be common to all plates. 

[43] A common succession of a fast (220—200 Ma, not 
documented in this paper) phase followed by a general 


Figure 7. Mean poles computed in successive 20 Myr 
time windows, every 10 Myr, using the common site 
longitude poles of Figure 6. Full radial scale is 10° only. 
Notice the small angular distance of each mean from the 
pole (on the order of 2° on average); 95% confidence circles 
are not shown, for clarity; all but two include the 
geographical pole. 
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a Europe 180° 


b North America 180° 


204 


90° 


Figure 8. (a) Europe and (b) North American synthetic 
APWPs (data in Table 4). Each pole corresponds to a mean 
pole computed by using a 20 Myr sliding window. Equal- 
area projection. Ages are in Ma; the values shown are the 
actual mean ages derived from the data in the appropriate 
window. 


standstill of Pangea at 200—180 Ma is followed by the 
breakup of Gondwana and the opening of the central 
Atlantic Ocean. An animated sequence derived from the 
APWP (www.ipgp.jussieu.fr/--besse/tpw jgr01) shows a 
sequence of general rotations around 160—140 Ma. The 
next phase of breakup of Gondwana results in the 110—100 
Ma events. The directional change of IND and AUS at 80 
Ma does not show as a major event. 


5. True Polar Wander 


[44] As was previously done by Livermore et al. [1984], 
Andrews [1985], ourselves (BC91), and recently Prévot et al. 
[2000], we have attempted to combine apparent polar 
wander based on paleomagnetic data (paleomagnetic refer- 
ence frame) with motions of hot spots with respect to 
lithospheric plates (hot spot reference frame) in order to 
derive true polar wander (TPW). Given the relative motions 
between the plates, the position of these plates with respect 
to the rotation axis (provided by paleomagnetic results under 
the GAD assumption) and the relative motions between the 
plates and a reference frame attached to the mantle (if such a 
frame can be found), TPW is defined as the motion of the 
mantle with respect to the rotation axis. The latter is defined 
by hot spots, which are assumed to form an array of fixed 
points that provides the mantle reference frame. 

[45] A number of the assumptions made to determine 
TPW may of course be found erroneous, such as the 
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assumption of hot spot fixity. This is a long debated subject 
[e.g., Molnar and Atwater, 1973; Molnar and Stock, 1987; 
Tarduno and Gee, 1995; Steinberger and O’Connell, 1998; 
Di Venere and Kent, 1999; Koppers et al., 2001]. In their 
analysis, Müller et al. [1993] find that motions between 
individual hot spots in the Atlantic and Indian Oceans are 
less than a few km/Myr, 1.е., close to an order of magnitude 
less than typical plate motion. These Indo-Atlantic hot spots 
therefore do form a relatively coherent frame for that vast 
part of the Earth’s lithosphere (and underlying mantle) 
ranging from North America to India. 

[46] The major improvement of the present analysis over 
that in BC91, where Morgan’s [1983] hot spot tracks on the 
African plate were used, is the use of the revised model of 
Müller et al. [1993], at least for the last 130 Myr. This 
model combines a number of advantages such as inclusion 
of hot spot tracks from AFR, NAM, SAM, IND, and AUS 
and the same updated plate motions by Royer et al. [1992] 
that we used in transferring the paleomagnetic data from 
one plate to the next in order to build the synthetic APWPs. 
Therefore, the kinematics and chronologies are compatible, 
potentially reducing significantly sources of errors that 
could have affected the BC91 analysis. Unfortunately, only 
two hot spot tracks are available (Tristan da Cunha and 
Meteor/New England) between 90 and 130 Ma, implying 
potentially larger uncertainties for that period. 


a South America 180° 


» Australia 


» Antarctica 


Figure 9. Same as Figure 8 for the (a) South American 
and (b) Indian, Australian, and Antarctic synthetic APWPs 
(see Table 4). 
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Figure 10. Comparison of two synthetic Indian APWPs. Triangles correspond to mean poles computed 
every 10 Myr using a 20 Myr sliding window (Table 4); dots correspond to mean poles computed every 5 
Myr using a 10 Myr sliding window, and their associated 95% confidence ellipses are shown (Table 5). 
The dashed ellipse is the confidence interval of the ill defined 104.6 Ma pole. Ages are in Ma. Equal-area 
projection. 


[47] The resulting TPW curve is shown in Figure 11, with 
95% confidence circles, for the best (10 Myr) time reso- 
lution available (Figure 11а) and for the 20 Myr time 
resolution (Figure 11b) (also see Table 6). Points are shown 
every 5 Myr (10 Myr) and every other point is statistically 
independent. Besse and Courtillot [1991] compared BC91 
TPW to two previous determinations by Livermore et al. 
[1984] and Andrews [1985], and we refer the reader to their 
Figure 3. In Figure 12, we compare BC91 (20 Myr 


resolution), the new ВСО1 curve of Figure 11a, and that 
obtained by Prévot et al. [2000]. 

[48] The first 11 points of BCOI in Figure 11а, corre- 
sponding to the period 0—55 Ma are all in the same 
quadrant, between 4 and 9? away from the present rotation 
pole, with 95% uncertainties ranging from 2 to 7°. They are 
not statistically distinct from each other and therefore could 
correspond to a standstill. A mean position can be calcu- 
lated from all (102) data points in that time window: it is 
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Figure 11. True polar wander paths (TPWP) deduced from 
the hot spot model of Müller et al. [1993] going from the 
Present back to 130 Ma, and that of Morgan [1983] going 
from 130 back to 200 Ma, with associated 9596 confidence 
ellipses (shaded light gray). (a) For 10 Myr sliding window. 
The dashed ellipse is the confidence interval of the 104.6 Ma 
pole; (b) for 20 Myr sliding window. Ages are in Ma and 
correspond to the actual mean age of the data in the 
corresponding window. Data from Table 6. 


found to lie at 85.1°N, 153.3°E (Ags = 1.5°), significantly 
displaced from the pole. The youngest pole at 5 Ma (mean 
age 3.0 + 2.4 Ma) is at 86.4°N, 166.8°E (Ags = 2.5°), also 
significantly displaced from the pole but identical to the 
55 Ma overall mean. It seems that TPW may have been 
negligible for an extended period of about 50 Myr, but 
accelerated a few Myr ago, with a velocity on the order of 
100 km/Myr. Actually, the TPW path zigzags around the 
pole with no clear track as early as 75 Ma. There is another 
clustering of data points from 79 to 104 Ma, then the pole 
jumps to a different location, which is well constrained by 
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the 116 Ma pole. The 95% confidence intervals of the 100 
and 110 Ma poles are too large to outline with certainty the 
details of the track which is followed. The track extends to 
142 Ma, when it turns back on itself. The path back to 
200 Ma is again rather jagged although there is a tendency 
to drift away from the current pole. When 20 Myr time 
windows are used (Figure 11b), the path becomes somewhat 
smoother and a succession of a standstill at 160—130 Ma, a 
circular track from 130 to 70 Ma, and a standstill at 50—10 
Ma follow each other. The true polar wander rate between 
130 and 70 Ma averages 30 km/Myr. The BC91 TPW curve 
was somewhat smoother, with a large track between 50 and 
120 Ma, in part because it was based on 20 Myr averages. 
However, differences between the two are of second order, 
given the uncertainties. TPW amplitude was smaller in 
BC91, with the same general features (tracks and sandstill) 
from Present back to 110 Ma, with maximum differences of 
order 3—4° in position between 70 and 90 Ma. On the other 
hand BC91 had a 120—150 Ma loop (or cluster) up to 8? 
away from the new ВСО1 path. 

[49] It is interesting to compare the new ВСО1 TPW 
curve with that recently derived by Prévot et al. [2000]. The 
main difference between the two studies lies in the fact that 


Table ба. True Polar Wander Path for 10 Myr Sliding Window" 


Age, Ma N ХОМ) ФЕ) Age K 
3.1 30 86.3 1722 2.6 105.2 
11.9 21 85.5 153.7 3.1 107.3 
14.8 24 84.9 151 32 85.1 
19.6 16 82.5 143.7 4.6 67.1 
26 14 85.4 143.9 53 58.5 
30.2 13 83.6 185.3 52 65.9 
33.7 12 83.5 172.7 4.6 90.4 
39 8 814 156.1 72 63 
46.4 12 84.3 136.7 54 70.1 
49.9 17 84.7 1342 41 76.5 
55 22 85.8 154.9 41 58.2 
60.7 24 86 191.2 43 49 
64 24 85.1 232 3.6 68.7 
68.3 15 85.6 245.7 4.9 63.4 
754 10 89.4 154.6 73 46.1 
79.1 9 83.1 300.3 6.9 57.6 
83.4 5 82.1 303.7 9.7 63.3 
91.4 8 85.8 303.7 3.5 247.6 
94.1 8 84.6 313 6.6 72.4 

100.2 6 81.6 342.4 11.1 38.4 
104.6 4 81.9 306 25.2 15.8 
110.8 7 75.5 301.7 83 56.6 
116.4 13 75 296.7 2.7 248.1 
120.1 13 72 286.6 32 172.3 
122.8 10 714 284.6 3.5 194.5 
130.9 4 69.1 292.4 7.6 148.9 
134.3 5 69.3 292.6 5.7 180 
139.3 5 63.7 286.2 8.1 89.5 
142 3 61.9 279.1 12.3 100.9 
153.9 4 74.1 310.7 6.6 196.1 
155.7 7 71.8 309.2 41 214.3 
158.9 5 69 296.2 Т7 100.5 
168 8 63.7 293.8 9.8 33.7 
170.2 13 66.4 299.5 6.6 40.5 
176.7 13 68.7 318.6 74 34.9 
1784 12 67.5 321.6 7.6 33.4 
1824 9 68.4 324.5 7.6 46.5 
191.3 13 62.4 305.4 57 53.6 
194.4 14 60.2 301.8 4.6 74.8 
198.9 8 63.5 296 5.6 100.4 


а Apparent polar wander path for global (Indo-Atlantic) hot spots for the 
last 200 Myr. 
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Table 6b. True Polar Wander Path for 20 Myr Sliding Window* 
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samples per site (versus our 6), the rest being similar to 
ours. Also they rejected data with potentially insufficient 


Age, Ma N ACN) o(°E) A K 6 ао i 

: » averaging of secular variation (K > 100). Their final data- 
3.1 30 86.3 172.2 2.6 105.2 н ; 

83 54 85.7 161 2 955 base contains 118 poles, i.e., less than half of ours. The 
18.9 38 85 149.5 27 764  TPW they obtain is in general similar to ours (Figure 12). 
29.5 23 84.7 151.3 3.8 664 However, their 95% uncertainties (not drawn in their figure) 
40 24 84.1 157 3.4 798 are very large because of the small number of data points in 
52.2 31 85 1517 3.4 60 : : ; : 

595 45 864 1963 28 sgg most time windows. Like us, they find a succession of 
673 34 864 229.4 32 60 episodes, with a long and faster track separating two stand- 
77.9 14 87.1 291 6 454 still periods. The most intriguing result they obtain is a 
90 13 83.7 308.6 49 73.66 phase of very fast wander culminating around 115 Ma, 
97.6 12 83.7 309.5 6.9 40.7 o DUNT 

when 20^ of polar wander take place at velocities in excess 

113.6 17 76.7 298 4.5 64.6 > А 

1191 20 736 291.8 2.6 164 Of 500 km/Myr, i.e., more than ten times faster than the 
126.4 14 70.5 287 2.9 183.1 mean rate we obtain on the track from 130 to 60 Ma. This 
136.8 7 66.2 286.3 6.3 93.1 result is based on two sets of individual poles dated around 
7 2 s pen Sb 2. 2 100 and 120 Ma, and more specifically to the recordings 
1734 21 673 308.7 58 314 around 115 Ma provided by South African kimberlites (114 
178.8 18 68.5 317.6 53 431 № 118 Ma). It so happens that our database provides a 
189.7 23 63.8 309.1 4.5 461 . sequence of rather well determined poles, particularly at 
196.7 19 61.4 298.5 4 718 — 90 Ma (Ags = 4.9°; N = 13 poles) and 119 Ma (A95 = 2.7°; 


"Apparent polar wander path for global (Indo-Atlantic) hot spots for the 
last 200 Myr. 
Prévot et al. use only paleomagnetic data from magmatic 
rocks and restrict themselves to basalts and andesites, or 
intrusive rocks with a similar composition. Their selection 
criteria were at least 10 sites (versus our 6) with at least 5 
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М = 20 poles). The angular difference is 10.0 + 5.5? and the 
corresponding polar wander rate in this 28 Myr period is on 
the order of 40 km/Myr. We can attempt the same compar- 
ison between the 119 Ma pole and a 114 Ma pole (Ags = 
4.2°): the angular difference is 3.2 + 6.5° and is not 
resolvable. We therefore find no support for the episode 
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Figure 12. Comparison between three TPW Paths (95% confidence intervals not shown for clarity): 
dots, ВСО1 (this paper); diamonds and dashed line, Prévot et al. [2000]; and triangles, BC91 [Besse and 
Courtillot, 1991]. Crosses are at 10° intervals (see also Figure 11 and Table 6b). 
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of super fast TPW suggested by Prévot et al. [2000]. As 
stated by Hargraves [1989, pp. 1851, 1852] in his analysis 
of the dating of the 114 and 118 Ma kimberlites in South 
Africa, “it should be pointed out that accurate radiometric 
dating of these bodies is difficult” and “inconsistencies 
between the sequence of apparent pole positions and the 
available radiometric ages may in part be due to uncer- 
tainties in the age.” The 114 Ma pole belongs to a 118— 
142 Ma petrochemical group, and the 114 Ma pole 
belongs to a distinct 83—114 Ma group. According to 
Hargraves [1989] “There is no clear pattern of migration 
with age" and “if it is argued that identical poles mean 
identical emplacement and magnetization age, then. . .the 
detailed validity of some of the radiometric ages must be 
suspect." The remarkable hypothesis of Prévot et al. is 
therefore hardly supported by such uncertain data. 
Actually, Prévot et al. [2000] average their data in four 
subsets (0—80, 80—110, 110—140, and 140—200 Ma) in 
order to reduce uncertainties. At that level, our two 
determinations are in good agreement, but evidence for 
any fine scale behavior or superfast episodes is lost. 


6. Discussion 


[50] The 10 and 20 Myr resolution TPW curves of Figure 
11 provide our current best estimates of true polar wander 
over the last 200 Myr. We confirm earlier findings [e.g., 
Andrews, 1985; Besse and Courtillot, 1991] that true polar 
wander appears to be episodic in nature, with periods of 
(quasi) standstill alternating with periods of faster TPW. The 
typical duration of these standstill periods in on the order ofa 
few tens of millions of years (50 Myr). Typical polar wander 
rates during fast tracks are on the order of 30 km/Myr. Also, 
because of all the uncertainties in models of hot spot 
kinematics prior to 130 Ma (and even possibly prior to 
90 Ma), we feel it is not safe to place too much weight on 
behavior prior to that time. The major event since then 15 
therefore the end of the 130 to 60 Ma period of relatively 
fast polar wander, with a standstill (1.е., no or little TPW) 
from 50 Ma (actually because of larger 95% confidence 
circles, possibly 80—50 Ma) to 10 Ma. However, uncertain, 
evidence for the fact that Earth emerged from that standstill 
to enter a new period of fast polar wander in a different 
direction 10 Myr ago (3 Ma at the higher resolution; see 
previous section) is particularly interesting. That period 
would then still be going on. Actually, the inferred rate 
and azimuth of this recent phase of accelerated TPW are 
compatible with the historical values based on direct meas- 
urements: between 1900 and 1990, the axis of rotation has 
been moving at a rate of 135 km/Myr toward eastern Canada 
(281°E; arrow on Figure 13 from Hulot et al. [1996]). 

[51] A legitimate concern regarding the above conclu- 
sions is due to the fact that the analysis 1s not truly global, in 
that it fails to encompass the Pacific plate. Petronotis and 
Gordon [1999] have compiled an APW path (nine poles 
from 125 to 26 Ma) for the Pacific plate, with four poles 
based on the analysis of skewness of ocean crust magnetic 
anomalies, three on seamount magnetic anomaly modeling 
[Sager and Pringle, 1988] and two unspecified. Using the 
Pacific plate versus hot spot kinematic model of Engebret- 
son et al. [1985], we have determined a corresponding 125 
to 26 Ma “Pacific hot spot only" TPW curve (Figure 13). In 
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Figure 13. Comparison between the “Indo-Atlantic” 
(open dots and dashed line; this paper) and “Pacific” 
(squares and solid line) TPW Paths. The Pacific TPW is 
based on the Pacific APWP of Petronotis and Gordon [1999] 
and the plate/hot spot kinematic model of Engebretson et al. 
[1985]. The arrow shows the mean polar motion between 
1900 and 1990 [after Hulot et al., 1996]. 


that sense, the TPW estimate derived in this paper could be 
termed an *Indo-Atlantic" TPW. The “Pacific” and *'Indo- 
Atlantic" TPW curves are compared in Figure 13. This 
comparison is interesting because the data sets they are 
based on are entirely different and independent. Despite 
some significant differences to which we return shortly, it 1s 
worth emphasizing that the two curves are similar in shape 
(tracks, amplitudes, azimuths), particularly the 300? longi- 
tude trending track from 130 to 70 Ma, though the two are 
offset (in the same general direction) by about 7°. More 
precisely, the confidence intervals intersect (though points 
are not in the intersection) near 125, 90, 60, 40 and 30 Ma. 
The main differences occur near 80 and 70 Ma: the 82 and 
65—72 poles derived from Petronotis and Gordon [1999], 
and our poles at 77 and 67 Ma (which have moderate 
uncertainties), are clearly distinct. 

[52] We emphasize that there are ongoing debates on the 
validity of the data used by Petronotis and Gordon [1999] 
to construct their Pacific APWP, and also on the question of 
fixity of Pacific hot spots with respect to each other and to 
the Indo-Atlantic hot spots. For instance, Tarduno and 
Cottrell [1997] have determined the paleolatitude (based 
on inclination-only data from cores) of the 81 Ma old 
Detroit Seamount, which is part of the Emperor chain, not 
far from its northern termination in the Kuril Trench. The 
paleolatitude (36.2 + 6.9/—7.2°) is distinct from that based 
on the 81 Ma pole of the Pacific APWP [Sager and Pringle, 
1988; Gordon, 1983] which is on the order of 20°N. 
Tarduno and Cottrell [1997] exclude the possibilities of 
inadequate sampling of secular variation, bias due to unre- 
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moved overprints or off-vertical drilling. They point out the 
uncertainties encountered when building an APWP for an 
oceanic plate, such as the Pacific, solely from inversions of 
magnetic surveys over seamounts and/or analysis of skew- 
ness of marine magnetic profiles. Petronotis and Gordon 
[1999] evaluated the quality of their own skewness data, 
which they rank from A (best) to D (worst). The uncertainty 
in the 73Ma mean pole based on A quality data is 3 to 4 
times larger than that of the mean based on all data (which 
they prefer). Yet the A quality mean pole is compatible, due 
to its large uncertainty, with the coeval 72 Ma pole of Sager 
and Pringle [1988]. It is therefore not clear that the con- 
clusion of Petronotis and Gordon [1999] (namely, that the 
two are significantly different, when all data from A to D 
quality are used) can be accepted. 

[53] Following several authors [e.g., Parker et al., 1987; 
Parker, 1991], Di Venere and Kent [1999] argue that the 
reliability of the Pacific paleopoles based on either model- 
ing of seamount magnetic anomalies or determination of 
skewness of marine magnetic anomalies should be consid- 
ered suspect. They recall that both are prone to numerous 
biases and could yield errors in excess of 10° in the position 
of the mean poles derived from them. 

[54] Let us now review briefly the question of hot spot 
fixity. Norton [1995] has suggested that the famous 43 Ma 
Hawaiian bend was actually a *nonevent," i.e., indicated a 
change in motion of the Hawaiian hot spot with respect to 
the mantle rather than a change in Pacific plate motion. 
Koppers et al. [2001] have tested the fixed hot spot hypoth- 
esis for Pacific seamount trails. They use seamount locations 
to first determine stage Euler poles, which they then test 
against observed *°Ar/*?Ar age progressions. The stages 
tested are 0—43, 43—80 and 80—100 Ma. Koppers et al. find 
that the 0—43 Ma Hawaiian and Foundation seamount trail 
pair is the only one compatible with the fixed hot spot 
hypothesis. The 0—43 Ma Louisville/Hawaii, 43-80 Ma 
Emperor/Louisville and 80—100 Ma Magellan/Musician 
trail pairs would all require relative motions of order (or in 
excess of) 10 km/Myr. The 43—80 Ma Emperor/Line pair 
shows particularly large discrepancies, requiring motions of 
at least 30 km/Myr. However, we may note that the notion 
of stage poles does require minimum knowledge on the ages 
of the seamount trails being fitted, and therefore the deter- 
mination of a geometrical stage pole is not strictly inde- 
pendent from age information. Koppers et al. also publish, 
but do not discuss, stage poles going back to 140 Ma. 
Altogether, four out of these six stage poles correspond to 
durations of only 10 to 20 Myr and the notion of a stage, 
with regard to dating uncertainties, becomes fuzzy. The 
authors emphasize that uncertainties in age progression may 
not entirely reflect true geological uncertainties, such as 
prolonged volcanic activity, rejuvenation by younger hot 
spots and uncertainties in location of hot spot (e.g., volcano 
hysteresis or hopping, control by lithospheric faults and 
structures). On the basis of careful analysis of 14 Pacific 
seamount tracks, using updated age determinations and 
bathymetry, Clouard and Bonneville [2001] show that these 
Pacific seamounts are created by different processes, most 
being short-lived and certainly not related to deep-mantle 
phenomena. Only the Hawaii and Louisville chains qualify 
as long-lived hot spots that can be robustly tested for fixity. 
We conclude that inter-hot spot motions of order 10 km/Myr 
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(but not necessarily much more) of (the major) Pacific hot 
spots are a distinct possibility. 

[55] Di Venere and Kent [1999] have addressed the 
problem of relative motion between the Pacific and Indo- 
Atlantic groups of hot spots. Though they demonstrate that 
some motion must have taken place between West and East 
Antarctica, based on geological observation and paleomag- 
netic data, they conclude that this motion cannot account for 
more than 20%, and possibly as little as 4% of the 14.5? 
offset between the observed and predicted positions of the 
65 Ma Suiko seamount on the Emperor continuation of the 
Hawaiian hot spot track. They also discuss the integrity of 
the Pacific plate and the role of missing plate boundaries 
and errors in kinematic plate circuits, and find that they play 
a small role. DiVenere and Kent conclude that most of the 
apparent motion between the two main groups of hot 
spots is real, with an average drift of about 20 km/Myr 
since 65 Ma. 

[56] Despite criticism on the quality of skewness data and 
seamount poles, which have strong bases, and increasing 
evidence that there can have been ~10 km/Myr motion 
between major individual hot spots or hot spot groups, we 
find the first-order agreement between the Petronotis and 
Gordon [1999] Pacific and our Indo-Atlantic TPW curves 
shown in Figure 13 quite remarkable. We therefore propose 
that our BCO1 curve can be considered as a good first-order 
estimate of global TPW and that the frame of reference 
based on the surface traces of the major hot spots deforms 
only slowly, slower than plates and plate boundaries move 
with respect to each other. In that sense, these hot spots can 
be used as a frame of reference for the underlying mantle, 
regardless of their dynamics and depth of origin. Hence 
TPW appears to be a truly global phenomenon, with tracks, 
cusps, standstills, and more generally amplitudes and azi- 
muths which are now reasonably well determined. Cottrell 
and Tarduno's [2000] proposal that TPW has not exceeded 
5? for the last 130 Myr is not vindicated. 

[57] We now turn to evidence for episodes of fast to 
superfast TPW, which have been proposed by a number of 
authors. With due caution and suggestions of alternate 
explanations, Petronotis and Gordon [1999] see possible 
fast polar wander at 80—70 Ma. Sager and Koppers [2000] 
reexamined the 130—40 Ma segment of the Pacific APWP 
based on poles derived from magnetic anomaly modeling of 
^ Ar/?? Ar dated seamounts (and only seamounts, with no 
skewness data) and concluded that a rapid TPW episode (300 
to 1100 km/Myr) occurred in 2 to 5 Myr at about 84 Ma. The 
back and forth motion on Sager and Koppers' Pacific APWP 
between 73 and 117 Ma, and even more so on their TPW 
curve, are uncomfortable features that could well be linked 
to problems in data significance and statistical robustness. 
The 84 Ma motion, between the remote yet quasi-coeval 
poles 84W and 84E, in entirely canceled by previous motion 
from 93 to 84(W) Ma, and subsequent motion from 84(E) to 
73 Ma (1.е., the 73 and 93 Ma points are not far from each 
other) Taken at face value, Sager and Koppers 84W and 
84E poles would lead to TPW velocities in excess of 500 
km/Myr, whereas the 93 and 73 Ma poles lead to values on 
the order of 50 km/Myr. Very recently, Cottrell and Tarduno 
[2000] have re-analyzed the data on which Sager and 
Koppers base their 84 Ma episode of rapid TPW. In addition 
to pointing out modeling or magnetization uncertainties in 
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the seamount data, they test the TPW episode against the 
reference paleomagnetic record from the Umbrian Appe- 


nines (Italy) and conclude that the episode must be an 
artifact of spurious seamount data. 
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[58] We noted above the suggestion by Prévot et al. 
[2000] of the existence of an episode of superfast polar 
wander (500 km/Myr) at 114—118 Ma but concluded it was 
likely an artifact. Note that Sager and Koppers [2000] find 
no evidence for the Prévot et al. [2000] event, and vice 
versa. Tarduno and Smirnov [2001] argue, in a comment on 
the Prévot et al. [2000] analysis, that reliable, selected 
paleomagnetic data from mid-Cretaceous (90—125 Ma) 
granites from North America disagree with either paleola- 
titude predictions based on the hot spot reference frame or 
with the superfast (stepwise) TPW event of Prévot et al. 
[2000]. As correctly identified by Camps et al. [2001] in 
their reply, Tarduno and Smirnov [2001] failed to indicate 
that the three predicted “‘paleolatitudes” for North America 
are not independent. In Figure 14, we have plotted the 
“latitudes” predicted over the last 130 Myr (1.е., the time 
with better ocean kinematics and hot spot data) for test 
points respectively in North America, South Africa and 
India. Three latitudes are shown: “APW?” is simply the 
paleomagnetic latitude predicted by our synthetic APWP, 
i.e., by the world paleomagnetic data set, when plate kine- 
matics are integrated. We have seen that this was compatible 
with the very data originating from each respective plate or 
continent. Under the assumption of a GAD field, this 1s 
simply the geographical latitude. “HS” is the “latitude” 
predicted in a reference frame where the (Indo-Atlantic) hot 
spots would have remained fixed with respect to the rotation 
axis (1.е., по TPW). *TPW" is the latitude predicted by the 
TPW curve, 1.е., motion of the point, as if it were attached to 
the hot spot reference frame, with respect to the rotation 
axis. The three latitudes (of which only the first, APW, 
corresponds to the true geographical latitude) are respec- 
tively linked to three rotations: APW is linked to Q4, the 
rotation of the plate with respect to the Earth's rotation axis, 
HS is linked to О», the rotation of the plate with respect to 
the hot spot reference frame, and TPW is linked to Оз, the 
rotation of the hot spot reference frame with respect to the 
Earth's rotation axis. The three motions are linked through 
the simple equation (1): Q1 = Q2* Q3. Therefore, any of the 
three latitudes can be derived from the two others, and they 
are by construction mutually consistent. In practice, we use 
APW and HS to deduce TPW. Figure 14a confirms that the 
three latitude estimates happen to be similar, and rather 


Figure 14. (opposite) “Relative” latitudes calculated for 
three test points: (a) North America presently at 34.5?N, 
75°W; (b) South Africa at 20°S, 20°E; and (c) India at 40?N, 
75°E. In each case, three "latitudes" are given from Present 
back to 130 Ma: APW is the paleolatitude predicted from the 
synthetic APWP derived for that plate in the present paper; 
HS is the paleolatitude in the (Indo-Atlantic) hot spot 
reference frame, 1.е., assuming no motion of the hot spots 
with respect to the rotation axis [Müller et al., 1993]; TPW is 
the “‘paleolatitude,”’ had the site remained fixed with respect 
to the hot spot frame of reference, while that frame moved 
with respect to the rotation axis, following our global TPW 
curve. The three estimates are linked since the first rotation 
(APW) is equal to the combination of the two others (HS and 
TPW; see text). 20 Myr averaging windows are used, and 
points are shown at their actual mean ages, with 26 
uncertainties derived from the 95% confidence intervals. 
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constant back to 80 Ma in North America. From 90 back to 
130 Ma, the three estimates diverge. The observed paleo- 
latitude (APW) remains constant, as has been long noted 
from NAM paleomagnetic data (and as is again emphasized 
by Zarduno and Smirnov [2001]). However, the HS and 
TPW estimates smoothly diverge in an opposite sense by 
roughly equal amounts. It so happens that they compensate 
each other in equation (1). Therefore, contrary to what is 
stated by Tarduno and Smirnov [2001], but in agreement 
with Camps et al. [2001], there is no disagreement between 
the three curves. We also find no evidence, as stated above, 
for the sudden jerk in TPW near 115 Ma, which would need 
to be compensated by an opposite, fortuitously coeval jerk 
in hot spot latitude. We note in passing that this agreement 
of the synthetic path with paleomagnetic data derived from 
granites somewhat alleviates worries on the quality of these 
data, as argued by Zarduno and Smirnov [2001]. 

[59] Figures 14b and 14c display similar results for points 
on the African and Indian plates, in order to show that the 
respective behaviors of the three latitude estimates of course 
depend on the distance and azimuth under which the APW 
and TPW paths are seen. In the Indian case for instance, 
TPW results in very little change in latitude, at least back to 
100 Ma, and the APW and HS latitudes are virtually 
identical. This was shown from paleomagnetic data of the 
Reunion hot spot trace on the Indian and African plates 
during ODP leg 115 [Schneider and Kent, 1990a; Van- 
damme and Courtillot, 1990]. For the 100—130 Ma period, 
the HS latitude of India becomes more northerly than the 
APW latitude, whereas the reverse holds for South Africa. 
Such curves can be calculated for any point on any plate 
used in our database. They allow quick comparison of 
actual paleomagnetic data, predictions from the master 
synthetic APW, and values predicted from the hot spot or 
TPW curves (reference frames). Figure 14 serves to illus- 
trate that the characteristics of these three curves depend 
very much on the location to which they apply. 

[60] In conclusion, none of the several suggested super- 
fast events is based on sufficiently robust sets of observa- 
tions. It remains reasonable to assume that many of these 
features correspond to erroneous individual data or other 
sources of error. Only the recent phase of TPW since 10 Ma 
prompts us to accept that TPW velocities on the order of 
100 km/Myr can be maintained over periods of millions of 
years, although we have no specific geodynamic explan- 
ation for this event (which some authors associate with 
deglaciation and rebound). 

[61] Whether even faster velocities over shorter time- 
scales actually occurred cannot as yet be considered as a 
strong constraint that should be modeled in numerical 
experiments. Studies of TPW and mantle dynamics have 
entered a new phase with the advent of flow models, where 
seismic tomography is used to infer 3-D maps of density 
heterogeneities that drive flow in the viscous mantle. For 
instance, Steinberger and O'Connell [1997] calculate the 
degree 2 nonhydrostatic component of the geoid and derive 
inertia perturbations on Earth over the last 60 Myr. Their 
results are in reasonable agreement with our earlier BC91 
TPW estimates, though slightly larger. However, their 
model has smooth and regular TPW changes, rather than 
the episodic structure interrupted by standstills which we 
find. Richards et al. [1999] calculate polar motion for 
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different Earth models. Interestingly, isoviscous mantle 
models predict TPW rates much larger than observed, and 
a significant viscosity increase in the lower mantle is 
required to stabilize the large scale pattern of convection 
and bring TPW rates closer to observed values. Richards 
et al. [1999] find occasional inertial interchanges of polar 
axes with a duration of 20 Myr, due to avalanching in the 
lower mantle (though only one such event occurs in a 600 
Myr numerical run). 

[62] None of these early models actually feature litho- 
spheric plates. The more recent study of M. Greff-Lefftz 
and P. Bunge (personal communication, 2000) and Greff- 
Leffiz [2001] explores the respective and cumulative effects 
of lower mantle viscosity, upper to lower mantle phase 
transitions and heat flux from the core on TPW estimates, 
combining 3-D spherical mantle circulation models with 
solutions to the equations of conservation of angular 
momentum. M. Greff-Lefftz and P. Bunge (personal com- 
munication, 2000) and Greff-Leffiz [2001] confirm that 
isoviscous mantle convection models predict TPW rates 
going from 100 km/Myr (the maximum acceptable value 
according to our study, applying to the current period) to 
1000 km/Myr and more (i.e., unacceptable values for us, as 
long as the elusive superfast episodes are not confirmed). 
Greff-Lefftz and Bunge find that phase transitions have 
little effect, and that combination of very high lower mantle 
viscosity (100 times that of the upper mantle, identical to 
the average they use for the lithosphere) with 12% bottom 
heating from the core [Davies, 1988] results in calculated 
TPW closest to that which we observe. “Because the 
rotation axis can only change as fast at the Earth's rotational 
bulge relaxes by means of viscous flow," inertial inter- 
change almost never happens. However, M. Greff-Lefftz 
and P. Bunge (personal communication, 2000) and Greff- 
Leffiz [2001] insist that their models strongly underestimate 
the true vigor of mantle convection, hence strongly under- 
estimate TPW rates. 

[63] Therefore, all current modeling still fails to some 
(sometimes large) extent to account for the slow values of 
typical TPW velocity (30—100 km/Myr), and even more so 
to account for the prolonged (~50 Ma) periods with almost 
no TPW (standstills). Also they predict rather smooth 
evolutions, rather than the alternating episodes which we 
feel we uncover from the data. The remarkable similarity 
between TPW estimates for the Pacific plate and the rest of 
the world, which are based on completely different and 
independent data sets, lends support to the idea that sig- 
nificant TPW, on the order of 10? or more, occurred since 
before the Cretaceous. The importance of the Pacific plate 
and severe limitations on presently available data from that 
plate point to the need for many more direct (paleomagnetic 
core) measurements as opposed to indirect/remote sensing 
determinations of magnetization direction (i.e., “skewness” 
or “seamount” data). The possible links between episodes, 
or major changes between TPW episodes, and either plate 
motion and plate boundary reorganizations or avalanches, 
plumes and other major geodynamical events occurring in 
the lower mantle should remain the topic of fascinating, 
ongoing studies. 
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